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SUMMARY 


i 


In  a  previous  report (>a>  s  statistical  theory  of  radar 
detection  was  presented  in  outline  form.  The  Mathematical 
details  were  omitted,  in  order  that  the  main  ideaa  and  results 
night  be  made  available  as  soon  as  possible. 

This  report  contains  the  mathematics  that  led  to  the 
results  presented  in  Ref. 28. 

In  addition,  several  subjects  are  briefly  discuaaed 
that  were  not  covered  in  Ref. 28.  Ihese  are  collapsing  loss, 
antenna  bema  shape  loss,  the  effect  of  signal  injection,  limiting 
loss,  end  moving  target  indication. 


For  ioforoacoa  a«a  pago  111. 
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SYMBOLS 


be  ob 


amplitude  of  sine  wave  relative  to  R.M.S.  noise  level 
one  of  the  independent  variables  in  the  function  Q( a, /3) 
ith  central  standard  moment 

one  of  the  independent  variables  in  the  function  Q(a,/3) 

half-power  antenna  beamwidth 

coefficient  in  the  Gram-Charlier  series 

characteristic  function 

delta  function 

base  of  natural  logarithms 

frequency 

confluent  hypergeometric  function 

Campbell  and  Foster  notation  for  characteristic  function 
Campbell  and  Foster  notation  for  anticharacteristic  function 
the  gamma  function 

probability  that  the  sum  of  N  noise  variates  will  exceed  the  bias  level 
i*h  Hermite  polynomial 
index,  subscript,  or  /-I 

incomplete  gamna  function  as  defined  by  Pearson  (8) 
modified  Bessel  function  of  the  first  kind 
Bessel  function  of  the  first  kind 
Ith  cumulant 

standard  Ith  cumulant,  or  sometimes  a  modified  Bessel  function  of  the 
second  kind 

integration  loss  # 

collapsing  loss 

generalized  Laguerre  polynomial 

number  of  excess  noise  variates  integrated  with  N  signal  plus  noise  variates 
false  alarm  number 


1  number  of  variates  integrated 
iyabol  has  a  different  meaning  in  HA- 15061. 

lyobol  ia  used  in  more  than  one  aenae  in  various  places,  but  other  meanings  should 


p* 

p 

p** 

p. 


<p‘ 

Q(a,0) 

ft” 

ft„ 


s 

<7* 

r* 
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«-(/) 

CD 

X 

y * 
K” 


to>  *  0)1/  -  1 

sine  wave  amplitude 
probability 

probability  that  noise  will  exceed  the  bias  level  at  least  once  within  false 
alarm  time 

probability  that  the  sun  of  N  variates  of  signal  plus  noise  will  exceed  the 
bias  level 

i th  derivative  of  the  error  function 
R.V.S.  noise  level 
modified  Lommel’s  function 
envelope  amplitude  or  radar  range  in  fi/fi0 
idealized  radar  range 

collapsing  ratio,  ratio  of  total  number  of  variates  integrated  to  those 
containing  signal 

cathode  ray  writing  speed 
standard  deviation 

(y*7)/cr>  semi -independent  variable  in  Gram-Charlier  series 

incomplete  Toronto  function 

i tb  moment  about  the  mean 

Lommel's  function 

normalized  envelope  amplitude 

i th  moment 

power  spectrum 

2  rrf 

power  signal - to-noise  ratio 
normalized  detector  output 

integrator  output  for  the  sum  of  N  variates 
bias  level 


*Thi«  symbol  has  a  different  meaning  in  RA- 15061. 

**This  symbol  is  uaed  in  more  than  one  sense  in  various  places,  but  other  meanings  should 
be  obvious. 


A  STATISTICAL  THEORY  OF  TARGET  DETECTION 
BY  PULSED  RADAR:  MATHEMATICAL  APPENDIX 


BASIC  FORMULAE  RELATING  TO  THERMAL  NOISE 


Both  th«  thermal  noise  voltage  across  a  resistor  and  the  noise  voltage  due  to 
the  shot  effect  in  a  vacuum  tube  approach  a  normal  distribution  when  the  number 
of  electrons  involved  per  second  in  the'  processes  tends  toward  infinity.  In  prac- 
tice,  it  may  usually  be  assumed  that  the  total  noise  voltage  between  any  two  points 
due  to  any  combination  of  thermal,  shot,  and  cosmic  noise  sources  can  be  represented 
by  the  distribution  function 


dP 


(1) 


where  4>Q  is  the  mean  square  value  of  the  noise  voltage* 18 This  distribution  is 
valid  provided  all  elements  involved  in  the  composition  of  the  total  noise  voltage 
have  been  linear. 


If  such  noise  is  now  passed  through  a  linear  filter  of  center  frequency  f , 
having  a  pass  band  which  is  narrow  compared  to  fm  ,  the  output  will  have  an  enve¬ 
lope,  which  has  a  probability  density  function 


dP  • 


(2) 


For  reference*  see  pege  111. 
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where  B  is  the  amplitude  of  the  envelope  and  \pQ  is  the  mean  square  noise  voltage, 
given  by 


u(f)  is  the  so-called  power  spectrum  of  the  filter  and  is  simply  the  square  of  the 
absolute  value  of  the  amplitude  transfer  function  of  the  filter. 

If  the  input  to  a  filter  consists  of  a  sine  wave  of  frequency  fm  ,  as  well  as 
noise,  then  the  probability  density  function  of  the  output  envelope  is* 


where  P  is  the  amplitude  that  the  sine  wave  would  have  at  the  output  of  the  filter 
in  the  absence  of  noise,  and  is  a  modified  Beasel  function  of  the  first  kind 
(see  footnote,  page  13), 

The  envelope  of  the  output  ha*  a  correlation  time  which  is  approximately  equal 
to  the  reciprocal  of  the  bandwidth  of  the  filter.  In  simple  language,  it  is  im- 
probable  that  the  envelope  will  change  by  an  appreciable  percentage  in  times  much 
less  than  the  correlation  time,  but  it  is  quite  probable  that  it  will  change  by  a 
goodly  percentage  in  times  large  compared  with  the  correlation  time.  It  is  probably 
a  good  approximation  to  assume  that  values  of  the  envelope  1/A f  seconds  apart  are 
independent,  where  A/  is  the  bandwidth  of  the  filter.  By  assuming  such  a  discrete 
process  it  is  possible  to  materially  simplify  calculations  which  would  be  very 
tedious  if  exact  integration  processes  were  used,  while  at  the  same  time  sufficient 
accuracy  is  obtained  for  most  practical  purposes. 

A  further  justification  for  this  assumption  in  the  pulsed  case  shows  in  the 
results.  Changing  the  factor  1/A/  to  fc/A /  for  the  correlation  time  has  only  the 
effect  of  changing  the  false  alarm  time  by  the  factor  k,  The  probability  of  detec* 
tion  turns  out  to  be  a  very  insensitive  function  of  the  false  alarm  time,  so  that 
if  k  is  any  factor  of  the  order  of  magnitude  of  unity,  the  results  are  affected  to 
a  negligible  extent. 


*  It  is  of  some  interest  to  note  that  the  seme  fora  of  dietribution  function  occurs  in  other 
probleaa.  For  inatence,  if  4ig  represents  the  mean  square  velocity  of  a  gas  dua  to  ordinary 
turbulence,  tod  P  represents  the  tranalational  velocity  of  the  whole  seas  of  gas  relative 
to  some  fixed  reference,  then  the  deneity  function  of  Eq.(4)  givee  the  probability  that 
the  total  vector  velocity  at  any  point  in  the  ges  will  have  a  magnitude  between  A  end 
A  v  dA(a4*. 

The  same  deneity  function  also  represents  the  probability  that  a  bomb  will  hit  at  a  distance 
between  A  and  A  ♦  dA  from  a  given  point  if  it  is  initially  siosd  at  a  point  whose  distance 
from  the  given  point  is  P.  The  mean  square  aiming  error  ie  represented  by  the  distri¬ 
bution  being  assumed  Gaussian  1471  . 
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DEFINITION  AND  EFFECT  OF  DETECTION 


A  detector  is  defined  as  any  device  whose  instantaneous  output  is  a  function 
of  the  envelope  of  the  input  wave  only.  Thus 


(-4) 


y  •  f\ 


where?  is  the  output  of  the  detector  and  v  is  the  normalized  amplitude  of  the  envelope. 
If  P/v'sZ'T  is  replaced  by  a,  Eq.  (4)  may  be  written 


dP  ■  ve  2  IQ{av)dv,  v>0 
dP  *  0  v <0 


Eq.  (5)  solved  for  v  is  v  ■  g(y) 


If  v  is  eliminated  from  (6)  and  (7),  an  equation  of  the  form 


dP  ■  f(a,y)dy 


is  obtained  for  the  probability  density  for  the  normalized  voltage  at  the  output 
of  the  detector  which  has  the  characteristics  given  by  Eq.(S).  For  example,  if  y 
•  vJ/2,  then  Eq.(8)  becomes 


dP  ■  e~r~*  I0  (2  Vxy^dy,  ?>0 


dP  •  0 


where  aJ/2  has  been  replaced  by  x,  The  quantity  x  may  be  identified  with  the  power 
signal-to-noise  ratio,  cortmonly  used  in  radar  literature. 


EFFECT  OF  VIDEO  AMPLIFIERS 

Since  a  complete  radio  receiver  usually  has  one  or  more  stages  of  video  ampli¬ 
fication  following  the  detector,  it  would  seem  that  one  would  want  to  calculate 
the  probability  density  function  for  signal-plus-noise  at  the  output  of  the  video 
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amplifier.  This  can  be  done  theoretically,  as  has  been  shown  in  an  excellent  paper 
by  Kac<as>,  but  the  mathematical  labor  is  great.  On  the  other  hand,  it  has  been 
shown  experimentally <aa)  that  the  signal  threshold  is  practically  unaffected  by 
the  video  bandwidth  until  it  becomes  less  than  about  lA  of  the  IF  bandwidth.  Since 
video  bandwidths  less  than  H  cf  the  IF  bandwidth  are  quite  uncommon  in  practice,  it 
appears  best  for  the  sake  of  simplicity  to  assume  the  video  bandwidth  infinite  in 
all  the  work  which  follows. 

When  the  results  have  been  computed,  assuming  an  infinite  bandwidth,  it  will 
be  possible  to  modify  them  in  an  approximate  manner  so  that  they  become  valid  for 
any  video  bandwidth.  This  is  explained  on  page  S9,  under  the  title  "Collapsing Loss”. 


PROBABILITY  OF  DETECTION  WITH  NO  INTEGRATION 

The  calculations  necessary  to  determine  the  probability  of  detection  when 
exactly  one  correlation  interval  is  available  are  quite  simple  compared  with  the 
case  where  the  output  over  many  correlation  intervals  is  available,  and  hence  the 
former  case  is  taken  up  first.  In  a  pulsed  system  this  corresponds  to  using  a  single 
pulse,  while  in  a  c-w  system  it  is  equivalent  to  observing  the  output  for  a  time 
t  m  \/bf ,  where  A/  is  the  over-all  effective  bandwidth.  In  either  case  this  amounts 
to  observing  the  receiver  output  for  one  correlation  interval.  If  the  output  exceeds 
the  bias  level,  the  signal  is  observed  or  detected  (see  pages  9-14  of  RA45061. 
A  Statistical  Theory  of  Target  Detection  by  Pulsed  ftadar<a#>,  hereafter  referred  to 
as  No.l,  for  complete  definitions  of  detection  and  bias  level). 

It  will  now  be  shown  that  the  probability  of  detecting  a  given  signal  x  is 
independent  of  the  detector  function,  everything  else  being  held  constant  and  only 
one  variate  being  taken  from  the  density  function  of  Eq.(8).  The  false  alarm  time 
has  been  defined  as  the  time  in  which  the  probability  is  %  that  the  noise  alone 
will  not  exceed  the  bias  level  (Eq.(15),  No.l),  but  it  will  be  best  here  to  keep 
things  general  and  denote  this  probability  as  Pg  ,  rather  than  as  %.  Eq.(15),  No.l, 
then  becomes 


P  »  1  —  pN/n  B  p 
A  0 


(10) 


where  the  subscript  N  denotes  the  number  of  variates  and  T  is  simply  an  abbrevia¬ 
tion.  From  Eq. (8) , 


S’  / 

^  Vi 


f<«) 


n  ■  I  /(0,y)  dy 


(11) 


where  the  symbol  y ^  is  now  used  for  the  bias  number.  Then  the  probability  of  detec¬ 
tion  is 


but  since  y  «  F(v),  or  v  «  g(y) ,  Eq .(11)  may  be  written 


and 


Therefore  Eq.(12)  becomes 


(13) 


(14) 


r 


ve 


IQ(av )  dv 


l0«.  ^ 


(15) 


which  is  independent  of  the  detector  function. 

The  integral  of  Eq. (15)  must  be  evaluated  by  approximate  methods.  This  function 
will  appear  in  several  places  subsequently,  and  is  defined  as* 


Q(a,/3) 


IAav )  dv  . 

0 


(16) 


Footnote  on  Q  Functions 

*  It  does  not  appear  possible  to  express  the  Q  function  in  terms  of  s  finite  number  of 
known  functions.  The  Q  function  is  similar  to  Lommel’s  functions  and  in  fact  can  be  ex- 

presaed  as 

,£±£lr  -l 

3  [iU^-i/3 *vae)-Ua{-i/3*,ia0)\ 

where  Ui  end  Uj  are  Lommel's  functions  of  the  first  kind.  This  identity  may  be  proven 
using  the  definite  integrals  given  in  WataonU),  pages  540  and  541,  especially  Eq. 5  of 
page  541.  By  successive  integration  by  parts,  the  Q  function  may  be  expanded  in  infinite 

series  giving 

Q(a,/3)  *  e 

r*®/  0\r 

Q(a.fi)  •  l  -  e  3  h  {a0)  • 

The  similarity  of  the  first  of  these  expansions  to  the  series  for  (/j(w,z)  given  in  Eq.(l), 
page  537  of  Watson,  is  interesting.  A  simple  expression  for  Q(a,a)  analogous  to  Eqs.(9) 
and  (10),  page  538  of  Watson  is 

Q(a,a)  ■  Tl  +  e"“8  J0(a3)l  (Continued  on  next  page.) 
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In  terms  of  this  notation,  Eq.(lS)  may  be  written 


p,  -  <?(«,/ 2  log, i/rt ) 


(17) 


This  is  the  probability  of  detection  only  if  Tj  »  l//p  where  is  the  time  available 
for  detection.  In  general  the  probability  of  detection  is  given  by  Eq.(18),  No.l, 


PH.y 


l-d-p,}- 


(IS) 


which  follows  from  the  definition  of  detection  given  on  page  9,  No.l. 
subscript  notation  PM  „  is  used  here  in  place  of  P'.  If  N  •  y,  P„  „ 

in"*'  " 

simply  as  . 


The  double 
is  written 


P,  r  (fl//?0)  can  be  calculated  by  means  of  Eqs.(18),  (17),  the  tables  of  Q,  and 
the  simple  relation 


1  2i/4 

*o  x 1/4  fa 


(19) 


(see  Eqs.(10),  (11)  and  (12),  No.l). 


Footnote  on  Q  Function*  ( Cant'd ) 


which  is  useful  in  special  case*.  An 
given  by  Rice(1,),  page  109: 

<?( a.y3)  »  i[i  -  <T‘(£-a)]  + 


asymptotic  expansion  for  Q  which  is  of 


(J-a)ar 

- 1 - «'  a  1 

2a  yjn 


fi-a  i+(/3-a)a 

♦a  +  *a*  ** 


value 


where  *"l(T)  is  given  by  the  error  function  of  Eq.(lOO).  This  expression  is  most  useful 
in  the  region  where  ay9  »  1  end  a  »  |/3-a| . 


The  Q  function  is  a  special  case  of  the  incomplete  Toronto  function  described  in  the 
footnote  on  psge  28  .  The  relstion  is 


<?(a,^) 


The  Q  function  is  grsphsd  in  Figs. 13  and  14. 

A  table  is  available  in  Ref. 47  but  the  intervals  are  too  Urge  to  be  of  general  use.  Project 
RAND  is  computing  an  extensive  table  of  the  Q  function  which  will  be  published  as  s  sepa¬ 
rate  report. 


GENERAL  CASE  -  INTEGRATION  OF  N  INDEPENDENT  VARIATES 

If  the  output  of  the  receiver  (or  filter)  can  be  observed  for  a  length  of  time 
much  greater  than  one  correlation  period,  it  ia  of  advantage  to  integrate  the  out* 
put.  The  simplest  concept  of  an  integrator  is  a  device  which  linearly  adds  the 
voltage  output  of  N  samples  from  the  detector.  The  time  elapsing  between  samplings 
must  be  at  least  one  correlation  period,  in  order  that  the  samples  may  be  considered 
to  be  independent.  If  the  sum  of  N  variates*  of  signal-plus*noise  exceeds  the  bias 
level  calculated  from  the  probability  density  function  for  N  variates  of  noise 
alone,  then  the  signal  is  said  to  be  detected. 

*  Readers  with  a one  atatiatical  experience  will  recognite  that  here  is  a  case  of  testing  a 
statistical  hypothesis.  It  ia  known  that  the  n  observations  yt ,  ,  -  -  -  -  y  cone  from 

a  universe  whose  density  function  f(y ,o)  depends  on  the  unknown  parameters  a;  it" is  required 
to  decide,  on  the  basic  of  these  observations,  which  of  the  two  values  a,  or  Og  ia  a  better 
estimate  for  a.  If  ia  the  true  value  of  a,  1st  pt  be  the  probability  of  asking  the 
mistake  of  choosing  aa  as  the  correct  value;  similarly,  if  aa  ia  the  true  value,  let  pa  be 
the  probebility  of  choosing  at.  Suppose  p.  ■  .05.  Than  a  atatiatical  decision  method  can 
be  devised  for  which  pt  ■  .05  and  for  which  p2  will  be  less  than  for  any  other  method  with 
the  tame  pt .  See,  for  example,  Kendall,  vol.2>  pp. 272*27 5 (®> . 


» 


The  integrator  may  take  the  sum  of  the  squares  of  the  N  variates,  or,  in  general, 
the  sum  of  N  variates  where  each  variate  has  been  processed  by  some  general  func- 
tion.  A a  long  as  the  same  weight  is  applied  to  each  variate,  the  integrator  will  be 
called  linear.  The  function  which  the  integrator  applies  to  each  variate  will  be 
called  the  law  of  the  integrator.  Any  nonlinear  integrator  will  be  inferior  in  opera* 
tion  to  a  linear  integrator  with  the  same  law  and  would  ordinarily  never  be  used 
intentionally  in  practice.  Cathode  ray  tubes  are  nonlinear,  however,  and  thus  fall 
short  of  other  types  of  linear  integrations. 

The  law  of  the  integrator  acts  in  exactly  the  same  way  as  the  law  of  the  detec* 
tor.  Thus,  if  the  detector  output  is  y  *  {F(v)  as  given  by  Eq.(5),  the  integrator 
output  is 


Y  -  ££(y)  -  £<t>lF(v)]. 

i  i 


(24) 


It  is  obvious,  as  far  as  the  theoretical  problem  is  concerned,  that  the  only 
function  of  importance  is 


« P(v)  «  <£[/'(»')]  . 


(25) 


There  will  be  an  infinite  number  ot  combinations  of  0  and  F  functions  which  will 
produce  the  same  function  0  and  hence  the  same  theoretical  results.  In  all  the  work 
that  follows,  the  output  of  the  combination  of  integrator  and  detector  for  one 
independent  variate  will  be  called  y  ■  0(v),  and  the  sum  of  N  variates  will  be 

.v 

Y  -  £y  .  (26) 

i 


The  symbolic  solution  for  the  case  of  N  variates  corresponding  to  Eq. (15)  for  one 
variate  is  not  too  difficult  to  obtain.  The  starting  point  is  Eq.v8)  for  the  proba¬ 
bility  density  function  for  one  variate.  The  characteristic  function  for  this 
distribution  is 


Xco 

/(  a,  y)  e*"*  dy  .  (27) 

OD 


The  characteristic  function  for  the  probability  density  function  for  the  sum  of  N 
independent  variates  is  then  simply 
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(28) 


■hWSu  1 
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j  ifc  *\j  ar-j  ar,  ■  wr->  m n-  >  w»,  * 


v  o  v-n  m  ka  v1  kakaka  um n  vi  *sn  v\ v  iu\ 


Corresponding  to  Eq.(ll)  is 
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and  to  Eq. (12). 


PN  •  f  G(a,JV 
Jh 


,Y)dY 


If  Y^  is  eliminated  from  Eqa.(31)  and  (32),  there  results  a  solution  for  ^  as  a 
function  of  ,  N,  and  a,  which  is  the  desired  result. 

It  is  found  in  most  cases  that  the  integrations  required  in  Eqs.(27)  to  (32) 
are  not  possible  in  terms  of  known  functions. 


THE  SQUARE  LAW  DETECTOR  WITH  /V  VARIATES 

It  seems,  by  a  process  of  trial  and  error,  that  the  best  possible  function 
for  t/r(v)  in  Eq.(2S)  to  produce  integrable  functions  in  Eas.(27)  to  (32)  is 


\ji(  v)  •  Av2  ■  y  . 


Though  this  represents  a  square  law  for  the  combined  detector  and  integrator  law,  it 
is  usual  to  think  of  it  as  representing  a  square  law  detector  coupled  with  a  linear 
law  integrator. 


I _ 


In  Eq.(33),  the  only  effect  of  the  constant  A  is  to  multiply  the  bias  level 
Y j  in  Eqs.(31)  and  (32)  by  A.  The  value  of  PH  in  Eq.(32)  is  independent  of  A,  It  is 
convenient  to  let  A  »  1/2,  or  y  ■  v3/2.  By  direct  substitution  from  Eqs.  (6)  and  (27), 


Xoo 

e'y~xI0{2  Sxy)e*y  dy  (34) 


where  *  •  oJ/2  and  p  ■  to>. 

This  integral  may  be  obtained  from  pair  655.1  of  Campbell  and  Foster*7*.  In  all 
pairs  taken  from  Campbell  and  Foster  it  is  necessary  to  replace  p  by  -p,  since  they 
use  eT*1*  for  the  first  integration.  As  long  as  the  same  notation  is  used  in  both 
directions,  the  order  of  signs  i*  immaterial.  In  order  to  avoid  confusion,  the  minus 
sign  will  be  used  in  the  exponent  in  the  first  transformation  and  the  plus  sign  in 
the  second  transformation.  Thus  all  of  the  characteristic  functions  which  appear 
hereafter  are  really  C(-p)  rather  than  C(p).  In  this  way  there  is  direct  agreement 
with  the  Campbell  and  Foster  tables  as  well  as  with  tables  of  the  Laplace  trans¬ 
form.  Equation  (34)  becomes 


C, 


(35) 


The  characteristic  function  for  the  sum  of  N  variates  is  then  simply 


e'*« 

(p+l)* 


(36) 


By  means  of  pair  650.0,  Campbell  and  Foster, 


the  probability  density  function  is 


",  ■  {£) 

-  o 


Y>  0 
Y<  0 


(37) 


Graphs  of  this  function  are  shown  in  Figs. 1-7.  The  density  function  for  noise  alone 
(x  ■  0)  is  found  most  easily  from  pair  431,  Campbell  and  Foster,  to  be 


BIAS  LEVEL  FOR  SQUARE  LAW  CASE 


The  bias  level,  Yb,  is  by  Eq.(31) 


r*Y“-le-r 

Jyk  <*-!>* 


dY 


(39) 


The  incomplete  gamma  function,  as  defined  by  Pearson <#>,  is 


I («*P) 


(40) 


In  terms  of  this  function,  Eq. (39)  becomes 


*lf  "  1  " 


(41) 


The  tables  of  I(/x,p)  extend  to  p  ■  50,  and  values  of  the  function  are  given  to  seven 
places.  Thus,  for  N  <  50,  and  n'  <  10* ,  the  bias  level  may  be  obtained  directly 
from  Pearson's  tables.  Other  methods  must  be  evolved  for  N  >  50  or  n'  >  10*.  The 
normal  approximation  to  Eq.(39)  is  unsatisfactory  for  N  less  than  several  thousand 
because  of  the  fact  that  the  integral  is  over  a  region  which  is  far  out  on  the  tail 
of  the  curve.  This  can  be  seen  from  the  Gram-Charlier  series  which  will  be  taken 
up  presently. 

The  integral  of  Eq.(39)  may  be  evaluated  directly  by  successive  integration 
by  parts  to  give 


r  .  £>-.r i.  sti .  fjtnjjtu  .... 

»  uv-u!  [  rl 


In  the  regions  of  interest  Y^>N»  1.  The  series  in  the  brackets  may  be  approximated  by 


1 


1 


(43) 
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in  this  region  so  that  Eq.(42)  becomes 


.,v-i  -K 


(N 


-1)!(lJr) 


NYU''' 

N'-(Yb’N+ 1) 


(44) 


By  the  use  of  Sterling's  approximation  for  iV!,  Eq.(44)  reduces  to 


jEni 

:y<>A i  *  lo®»  a?). 

U-rr 

(45) 


Though  the  expression  looks  more  cumbersome  in  this  form,  it  is  actually  much  simpler 
to  use  in  calculations  than  is  Eq.(44).  Substituting  for  from  Eq.(21)  gives  the 
expression 


log,0n  -  0.24  ♦  |  logw  N  ♦  log  10  ( yfc-iV+l) 
♦  0.434  ( 1T6-/V)  -  N  log10  $  . 


(46) 


Graphs  of  this  function  are  shown  in  Figs. 8  and  9.  For  N  •  1,  the  exact  expression 
for  Yb  from  Eq.  (39)  is 


Yb  •  2.3  logJ0  n  ♦  0.  37 


(47) 


whereas  Eq.(46)  for  N  ■  1  reduces  to 


Yb  -  2.3  logjQ  n  ♦  0.45 


(48) 


The  difference  is  seen  to  be  practically  negligible. 
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CUMULATIVE  DISTRIBUTION  FOR  N  VARIATES  OF 
SIGNAL  PLUS  NOISE  -  SQUARE  LAW  DETECTOR 

Knowing  the  required  bias  level  for  a  given  false  alarm  interval,  it  is  now 
necessary  to  integrate  the  density  function  of  Eq.(37)  from  this  value  to  infinity 
to  give  the  probability  of  detection  for  a  signal  of  strength  *,  thus 

4=1 

"(J.)  3  e-Y-NXIN_l(2^?)dY  .  (49) 


This  integral  is  not  soluble  directly  in  terms  of  well-known  functions.  The  order 
of  the  Bessel  function*  can  be  reduced  in  steps  of  1  by  successive  integrations  by 
parts,  so  that  the  last  remaining  integral  is  of  the  type  given  by  the  Q  function 


*  The  following  art  soma  of  the  useful  identities  concerning  the  modified  Bessel  functions  of  the 
first  kind:  + 

r.«  (if  P  -s  .s  1 


/,<*)  »  {-if  J*{  is)  ■ 


it) _ .  liL  I  „  _.«■  „  .i*_ 

r  .0  r!(ser)[  N 1  L  t-iiln+2H 


J0< s)  •  1  ♦  4j-  ♦  •— 

Asymptotic  expansion: 

/  <*>  I"1  * 

"  vTmf  L  1U8*)  J 

(?)" C  *>•■# 

«/'<■)  -  *!"{,)  *  s/-+i(s)  -  -«/„<*>  ♦ 

/„(  i)  ■  s) 

/  s"/. _,<*)ds  •  s"/„ (s) 

J"  ilg(  t)di  ■  i/t(  s) 

¥  v*>  -  w*>  -  w> 

/  e‘/0(x)d*  -  *e'Cf0(*>  -  Jjlx)] 

/  e“*f0(*)dx  •  *e"*[/0(s)  ♦  <*)] 

f  e*Ii(M)dM  ■  e*Cu-*)/0(s)  ♦  s/j  <*)] 

f  e"*/4(*)d*  ■  e'* [(!+*)/„(*)  ♦  s/^s)] 


Relations  between  the  fH  (unctions  and  the  hypergeoeatric  functions  will  be  found  in  the  foot¬ 
note  on  p.  21. 
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of  Eq.(?,6).  An  easier  way  to  arrive  at  the  same  result  is  by  the  use  of  the  charac¬ 
teristic  function.  To  get  the  cumulative  distribution  from  -00  to  Y  of  any  density 
function,  it  is  only  necessary  to  find  the  anticharacteristic  function  of  C/p,  where 
C  is  the  characteristic  function  of  the  given  density  function  (see  pair  210,  Campbell 
and  Foster).  Thus  from  Eq.(36), 


The  term  l/p(p*l)#  may  be  expanded  in  a  series 


p(p*l)N  p(p+l)  (p+l)a  (p+1)3  (p+1)* 


The  mate  of  the  first  termofthe  series,  by  pairs  210,  and  655.1,  Campbell  and  Foster,  is 


'  e-y I0(2/f?Zy)dy 


The  first  two  terms  of  f*  are  thus 


1  -J  a"y  -**!„  ( ly/Rxy)  dy 


/•« 

I  t/e'  3 

J/IrJ 


IAv/TNx)  dv 


Q{m~x,ASTh) 


using  the  definition  of  Q  from  Eq . ( 16 ) .  All  the  succeeding  terms  may  be  obtained 
by  using  pair  650. 0  Campbell  and  Foster, 


Nx  r-  l 

.p+I  /  v  v  “T* 


M*te  of  (7TT7  ia(£)  e‘yW2^) 


•  Aa  in  Caepbell  and  Foater,  /  ia  hare  used  in  place  of  u/2v  or  p/2w» . 
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From  Eqs.(53)  and  (54), 


Pv  -  Q(/Wxt/2Y^  e‘VW|5(^)a  IrJ 2/^)  .  (55) 


This  form  of  solution  for  ^  is  practical  for  numerical  calculation  only  where  N  is 
less  than  about  10. 

The  characteristic  function  in  Bq.(50)  can  be  expanded  in  another  manner  using 


l  .  y  1 

P(P  +  l)N  rr^y  (p*  Dr 


This  leads  to  an  expression  for  /-^complementary  to  that  of  Eq.(55)  of  the  form 

*.•  >-  «'v"  f  (&)*  t-,  •  <»> 

raN+l 


By  equating  (55)  and  (57),  one  obtains  one  of  the  known  expansions  for  Q  given  in 
the  footnote  on  page  5.  Equations  (55)  and  (57)  may  also  be  obtained  directly  from 
Eq. (49)  by  repeated  integration  by  parts.  Equation  (57)  may  also  be  obtained  di- 
rectly  from  Eq. (55)  by  means  of  the  identity 

tl  .  £  tnIn(x)  (58) 

h«-od 

given  in  McRobert<a>(  page  32,  and  one  of  the  known  series  for  Q> 

For  the  special  case  ^  ■  Nx,  the  function  Q  of  Eq.(55)  is  simply 

■  2  [l  *  e*J>,6lo(2ifr)]  (59) 

(see  footnote,  page  5  ),  and  Eq.(55)  becomes 

4  ■  J  ♦  *'"‘[^.(2'.)  *  1,(2! i)  *  «21i)  W21i)]  • 
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This  formula  is  useful  for  checking  special  points  for  values  of  N  around  10  or  below. 

None  of  the  methods  developed  above  are  suitable  for  calculating  PH  for  large 
values  of  N. 

In  the  next  section  the  general  method  of  Gram-Charlier  series  is  developed, 
which  will  be  useful  in  a  number  of  succeeding  problems  concerning  distribution 
functions  over  large  ranges  of  variation  of  N. 


EXPANSION  OF  FUNCTIONS  IN  GRAM -CHARLIER  SERIES 


The  function  0(y)  is  defined  by 


jr 


The  Hermite  polynomials  may  be  defined  by  the  relation 


Vs 

(-1)1  „2L 

'  bk*  '“i 


where  the  superscript  i  stands  for  the  i111  derivative  with  respect  to  y.  The  4> 
functions  and  the  Hermite  polynomials  are  biorthogonal ,  that  Is 


/*« 

H^yrfiy)  -  Si;.  -  0 ,i/j 

%/-  CD 


"  1 1  i  ■ 


Therefore  it  is  possible  to  expand  any  reasonable  function  in  a  series  of  the  form(8> 


f(y)  ■  ^  ^  (y) 


The  coefficients  a.  may  be  evaluated  in  a  manner  analogous  to  the  Fourier  series 
methods  by  multiplying  both  sides  of  Eq.(64)by  /f(y)  and  integrating  from  -®  to  a.  All 
terms  drop  out  but  one,  giving 


(y )f(y)dy  . 


S'1  * 


r  •  V*  iVXV'iy'/H 


*•5 

3 

S: 


l 

£3 


It  is  usual  to  make  the  substitution  t  *  y  -  y/cr  before  making  the  expansion,  thus 
causing  the  second  and  third  terms  of  the  series  to  vanish.  The  notation  is 

y  »  Vj  ■  the  average  value  of  y,  or  the  first  moment 
-  v*  ■  the  variance 

v  •  the  nt&  moment  of  the  distribution. 

ft 


ynf(y)dy 


Equation  (64)  is  replaced  by 


f(y )  •  g(  t )  -  y  e.01  (t ) 

i*o 


and  Eq.(65)  by 


ir-J 


It  follows  at  once  from  Eq.(68)  that  c„  ■  1 /cr,  »  c2  ■  0.  The  moments  about  the 

mean,  or  the  central  moments,  are  defined  by 


r 

«/.  go 


( y-y)fiy)dy 


and  the  standard  moments  about  the  mean  by 


a .  ■  — 

*  or* 


a  ru  m  >rs  xy  v*  •'■u  vuwvwn.  w  w w  w  ww^ry^y^Aj'vvjv^/wyv  awskiv*  aaa*  r*.  iwy*  />  r>  tv*  .r*nj*>v.K>.-Jw 


The  coefficients  c 
few  are 


Formulae  for  the  m' 


j  in  Eq. (68)  may  be  easily  written  in  terms  of  the  a’s.  The  first 


c4  •  7j(a«"3)  (71b) 

cs  “  (as“10aj)  (71c) 

c«  "  -7(a4"15a4+30)  (71d) 

c7  Z  "4  (ar-21a,+  10Sas)  (71e) 

*  • 

c«  “  7r(as“28a«+2l0a4‘315)  (71f) 

C»  "  (a»"36a7+3T8CS"1260a3)  *  (71g) 


s  in  terms  of  the  V  a  can  be  obtained  directly  from  Eq. (69),  giving 


(72a) 


M,  -  v3  ~  iv2vl  ♦ 


(72b) 
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M4  *  ^4  “  *vjvi  *  6vi v\  -  3l/i4 


(72c) 


Sv4Vj  +  lOi/jV,1  -  10v2vt3  + 


(72d) 


^5 


n6  ■  v6  -  6v$  Vj  +  15v4vf  -  20VjV3  +  lSVjv}  -  5^® 


(72e) 


Continuations  of  this  series  are  obvious. 

The  process  of  obtaining  the  Gram-Charlier  expansion  is  now  evident: 

1.  Find  the  moments  of  the  distribution. 

2.  Obtain  the  central  moments  from  Eq.(72). 

3.  Obtain  the  standard  central  moments  from  Eq.(70). 

4.  Obtain  the  coefficients  from  Eq.(71). 

5.  Write  the  series  for  /(>)  from  Eq.(67). 

It  turns  out  that  the  best  grouping  for  the  terms  of  the  series  of  Eq.(67)  is 
different  from  the  natural  sequence Such  a  regrouped  series  is  termed  an  "Edgeworth 
series"  and  is  actually  used  in  this  work.  The  grouping  used  by  Edgeworth  is 


0 

(73a) 

0, 

3 

(73b) 

0, 

3,  4,  6 

(73c) 

0, 

3,  4,  6,  5,  7,  9  . 

(73d) 

This  means  that  if  the  0  and  3  terms  are  used  as  the  first  approximation,  the  addi¬ 
tion  of  terms  4  and  6  gives  the  next  order  approximation,  and  so  forth. 


MOMENTS  OF  SIGNAL  PLUS  NOISE,  SQUARE  LAW  DETECTOR 

The  moments  of  a  distribution  may  be  obtained  by  using  the  characteristic 
function  as  a  moment  generating  function*.  Thus 


vi 


(74) 


•  Kend.il,  p.54(#) . 
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iwl  \  \, 


Jjtv. 


In  the  case  of  the  distribution  function  for  the.  sum  of  N  variates  of  signal  plus 
noise  with  a  square  law  detector,  the  characteristic  function  is  given  by  Eq . ( 36 ), 
and  the  moments  are 


H 


(-1)‘ 


e _ e 


dp1  (p+l)* 


(75) 


p»  o 


Though  the  first  few  moments  may  be  obtained  by  direct  differentiation,  it  is  better 
in  this  case  to  expand  in  a  McLaurin's  series  and  obtain  the  coefficient  of  p l/ i ! .  Thus 


JJL 


Sx 


(Nx)* 


(p+1)*  (p  +  1)"  (p  +  l>*+1  (p-1  -  2! 


The  coefficient  of  px/i\  is,  by  direct  expansion  of  each  term  in  Eq.(76), 


>ii£ 


‘•“  TFIirL1  *  *  “Turn) - IT*- 


■  (JV+ 1-1)1 

(‘1)l  OKI)!  «/»(JV+i,Ar,iVx) 


(76) 


(77) 


(78) 


where  jf,  is  the  confluent  hypergeometric  function*  Thus  the  moments  are 


*  The  following  are  some  of  the  useful  relationi  concerning  the  confluent  hypergeometric  function: 

F  (a,  c,  x)  -  ni£i  T  IfrttEl  -  1  ♦  fie  .  f,  («.  b.  e.f ) 

11  r<e)  r«o  r!  F(  c+r)  *  l!  e(c+u  2!  •  a  l'  o/ 


Asymptotic  expansion: 

.F  ( a,  e,  i)  *  XtdlL  fl  ♦  il=gli£=a2  ....] 

1  1  P(a)*e-*  L  s  1 !  J 

Runner's  first  transformation: 

jfj  ( a,  c,  t)  ■  e'jfj  ( c-a,  c, -*) 

Rummer's  second  transformation; 

*>  -«Vi("*i£) 

Recursion  relations: 

a^  (a+1,  e,  i)  ♦  ( a-e>  ^  ( a-1,  c,  i)  ■  ( a+2a-c)  f  ( a,  c,  i> 
( «>1,  c,  *)  ♦  ( e-a) s ^  (  a,  cel,  j)  ■  c(a+s)  (a,  c.  i) 
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Lim  6  “*  a> 


(Continmd  on  ntxt  pat*-) 


i  wn  is  a  ks.a  m_as_>a  aa  as  .sere  as  M  m  m-Un  Ll  s^ahwusasu  m"w  ar  w  ww  ■  vi  e  \j  eu  sueu  iar-u  w  V  ad  W  t#U  iS\y  W^>ar\j  trtJ  mr\j_jwu  WVi  nru  arfii  «r\i  awv  _ 


(79) 
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(/v+t“l) !  u .  . . 

*V  “  — )■;  e*"*IF1(/V  +  i.iV.JV*) 


UV+i-1)!  ,  .  „  „  % 


Hypargaomatric  function  ( Cont'd ) 

s^fa+l.c.  r)  ♦  ( 1-e)  (  a,  e-1,  i)  ■  <a41-e)  ^  (a,c,  i> 

•  /  i)  ♦  Jt^(a,c+l,j)  ■  Cj/f  (a.c.i) 

(a«e) /j  («”1,  e,  i)  ♦  <  c-1)  f  (  a,  e-1,  j)  ■  (  a+a-1)  ^  (a,  e,  *) 

(  c-a)  (  a,  e+1,  i)  ♦  e(  c-l)  Ft  (  a,  e-\,  z)  «  e( *+c-l )  jfj  (a,  c,  a) 

,fj  (a, c, i)  »  *  (a+1, e+l,  i) 

Relations  between  hypargaomatric  function*  and  other  functions: 
t/j  (a,  a,  a)  ■  «* 

^  (a,  «+l, -a)  ■  si“*  f0  a-1^ 

using  Psaraon'a  notation  for  ths  incomplete  gamma  function. 


^(1,0+l.j)  «  «-VT\.*l)I^,  a-l) 

/i(f  f’*“)  ■  Ko*  m$  "f  * 


tJJ  (-n,  1,  j)  ■  Lb(  *! 


(original  Laguarra  polynomial) 


.f.  (-n.a+l.i)  £(,) 

1  1  P  (a+l+n)  " 


(generalized  Laguarra  polynomial) 


.1  ()•»-)■  •'.(*)] 
/.H 

Mr  *•-•)■  ‘  *  -'.(f)] 


*r\nnmn  *m\  r*  \  «*•  >  fu  *-r»  u»n-*a*» 


Eq.(80)  being  obtained  from  Eq.(79)  by  Rummer's  first  transformation.  The  first 


four  moments  are 

v,  -  iV{  1+x)  (81a) 

vt  -  (iV*)J  ♦  2Nx (N+l)  ♦  N(N+ 1)  (81b) 

v3  •  (Nx)3  ♦  3(jV*)a((V+2)  +  3/Vx(/V+l)(JV+2>  ♦  N(N*l)(N*2)  (81c) 

v4  -  (Nx)*  +  *(Nx)3{N+ 3)  ♦  6(Nx)a  (A(+2)(iV+3)  (81d) 

♦  4Nx(N+l)  (N+2)  (iV+3)  ♦  JV(N*1)  (jV+2)  (N+3)  . 


The  generalised  Laguerre  polynomial  LjfH*)  is  defined  by 


.  .  Ha+l+n)  _  , 

1”(,)  *  TiruTT)  • 


(82) 


Comparing  (80)  and  (82),  it  is  seen  that  the  moments  expressed  in  terms  of  the 
Laguerre  polynomials  are 


v.  ■  ilL[s’lH~Nx)  , 


(83) 


Another  generating  function  for  these  polynomials  is  available  through  the  relation 


J1 

dtn 


(84) 


The  moments  about  the  mean  may  be  expressed  in  terms  of  the  moments  about  the  origin 
by  means  of  Eq«.(72a-e),  resulting  in: 


M0  "  1 


(85a) 


M|  •  0 


(85b) 


jj. 2  ■  2 Nx  +  N  ■  /V(2x+1)  ■  c r3  (85c) 

«  6iV*  *  2N  •  2N(3x+l)  (85d) 

M4  ■  12 ( /V* ) 3  +  12JVx(tf+2)  ♦  3W/V+2)  (85e) 

M,  ■  120  (Me) 3  ♦  20iVx(5IV+6)  *  4/V(5tf*6)  .  ( 85 f) 


A  generating  function  for  the  central  moments  may  be  obtained  by  multiplying  the 
generating  function  of  Eq.(75)  by  e,w» giving  (see  pair  207'  Campbell  and  Foater) 


(-I)*®"** 


sw- 


(86) 


The  momenta  of  Eqi.(85a-f)  are  moat  easily  obtained  by  logarithmic  differentiation 
in  Eq. (86) . 

The  standard  momenta  about  the  mean  are  obtained  from  Eq.(70),  and  are 


S  "  1 


(87a) 


a,  «  0 


(87b) 


a3  "  1 


(87c) 


2(3**1) 

l»  ’  1 
Nl/i(  2x+l) 3 


(87d) 


6  (  4x+ 1 ) 


a,  ■  3  +  ,  • 

4  JV(  2*+l)3 


( 87e) 
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There  is  an  approximate  method  of  computing  the  significant  part  of  a(  which  is 
based  on  the  fact  that  c#  of  Eq . ( 7 Id )  is  always  nearly  equal  to  c*/ 2  (see  page 
259,  Fry<*>).  Thus 


as  15  u4  ♦  10  a*  -  30 


(88) 


or 


a4  as  is  ♦ 


10{l08s*»78**13 ) 
N(  2*4.1)* 


(87f) 


For  noise  alone,  the  moments  are  given  by 


(A4i-1) ! 

Vi  ’  UV-1)! 


(89) 


and  the  central  momenta  by 


Mi 


(ft*i-l)! 

<(v-n! 


/,(-(,  l-i-w.-yv) 


(90) 


Equation  (90)  was  obtained  from  Eq.(86)  by  putting  *»0  and  expanding  in  a  series. 


THE  GRAM-CHARLIER  SERIES  FOR  THE  SQUARE  LAW  CASE 

The  coefficients  of  the  series  may  be  obtained  by  use  of  Eqs. (71a-d)  since 
the  standard  central  momenta  are  now  known  (Eqa.B7a-f).  They  are: 


eo  “  1 


(91a) 


C,  m  C3  m  o 


(91b) 


ej  "  ' 


(91c) 
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4  4iV(  2x  +  l)3 


m  ( ***  1 )  3 
C«  “  18N(2x  +  l  )*  * 


(91e) 


From  Eq.(67)  Che  required  aeries  is 


dP 


d>°(t) 


*  c40J(t) 


♦  c4<P*(t)  +  C#0*(t 


(92) 


where 


t 


(93) 


Vj  ■  JV( l+x) 


(94) 


(7  ■  /JV(  1+2*1) 


end  the  c's  are  given  by  Eq*.(91a-e). 

Note  that  the  grouping  of  terms  is  tccording  to  the  Edgeworth  scheme  given  in 
Eq.(73).  Note  further  that  as  N  tends  to  infinity,  all  the  coefficients  go  to  zero 
except  c0 ,  Thus 

.  .jziHi! 

dP  •  -=r-e  icr*  dy  (95) 

v2ft  c 


as  N  -  oo,  Jn  terms  of  iY  and  x 


Eq  .(95)  is  precisely  a  statement  of  the  central  limit  theorem,  and  the  derivation 
given  ia  essential ly  a  loose  proof  of  the  theorem. 

The  cumulative  distribution  is  easily  obtained  from  Eq.(92)  by  means  of  the 
simple  relation 


/■ 


(97) 


giving 


'  PN 


f(y)dy 


-c td>3(T) 


e4  <p*  (T)  -  c<0,(T)  ----(98) 


where 


and 


T 


(99) 


•“*  da 


(100) 


The  function  d>“1 (T)  is  tabulated  in  the  W.P.A.  tables This  differs  from  the 
definition  given  for  4>~l(y)  in  Fry,  page  456  but  is  used  here  because  of  the  W.P.A. 
tables. 

The  aeries  of  Eq . ( 98 )  was  used  to  calculate  all  the  curves  of  Figs. 1-50,  No.  1, 
with  the  exception  of  the  cases  where  N  ■  1.  In  most  cases  the  first  two  terms  of 
the  series  are  sufficient,  though  in  some  regions  of  small  P  four  terms  are  needed. 


SAMPLE  CALCULATION 

Assume  N  •  10,  n  •  10* 

From  Fig.  8,  or  Eq.41,  •  30.0 

B 

Let  /T  *  1.0  so  that  *  ■  1.0 
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From  Eq. (94) , 


v,«  NU+x)  m  20.0 
er  •  vWi* 2x)  ■  5.48 


From  Eq.(99), 


K-u  30.0  -20.0 

■  -V'— r 


1.830 


0-‘(T)  -  0‘‘ (1.828)  ■  0.9325 
j(l-0-‘(T))  -  0.0338 


From  Eq. (91c), 


&UL 


3^  (:*♦!)• 


-0.081 


From  Eq. (91d), 


'  °-om 


From  Eq. (91e ), 


cl  " 


.  0.0033 


l£W(2x+l) 


0s (1.828)  -  0.174 
4>J(  1.828)  ■  -0.470 
0*( 1.828)  •  0.990 
<:,<*>*( T)  -  -0.0141 
c4*»(T)  -  -0.0007 
c40*(T)  ■  +0.0032 

P  -  0.0338  ♦  0.0141  ♦  0.0007  -  0.0032  ■  0.0452 
This  point,  P  ■  0.045,  ■  1,  may  ba  found  on  Fig. 20,  No.l. 


(See  p.64  for  raferancea  on 
tables  of  tha  derivatives 
of  the  error  function.) 


INTEGRATION  LOSS,  SQUARE  LAW  DETECTOR 

It  ia  of  interest  to  express  the  effect  of  noncoherent  integration  as  a  loss 
with  respect  to  coherent  integration <41> .  This  may  be  done  by  defining  the  integration 

loss  as 

Nx. 

L(  ■  10  log, -  (100a) 
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wwvs m  a  d  s 


■  n  e  m  «  ■  w  w  «a  m  <a  a  a  »  w  ■  n  » a m wan mja  pan 


where 


JV  ■  number  of  pulses  integrsted 

*j  “  required  value  of  signal -to-noiae  ratio  to  produce  given  probability  of 
detection  for  N  »  1. 

»  required  value  of  signal-to-noise  ratio  to  produce  the  same  probability 
of  detection  for  N  ■  /V. 

Thus  is  a  function  of  P  and  n.  However,  it  turns  out  that  the  dependence  on  P 
and  n  is  very  small. 

In  the  case  of  coherent  integration,  *2  is  always  equal  to  x^/N,  so  that  Li  ■ 
0.  With  noncoherent  integration,  *r  is  always  greater  than  x./N,  so  that  noncoherent 
integration  is  never  as  efficient  as  coherent  integration.  The  results  of  calcula¬ 
tions  are  given  in  Figs. 10  and  11.  One  observes  that  the  dependence  of  Li  on  P  and 
n  is  quite  small.  Thus  by  means  of  the  graph  of  Fig. 12,  which  gives  *  as  a  function 
of  P  and  n  for  N  ■  1,  and  any  one  of  the  curves  of  Fig. 10,  it  it  possible  to  obtain 
a  fairly  accurate  value  of  x  for  any  P,  n  and  N. 


GENERAL  CURVES  OF  THE  CUMULATIVE  DISTRIBUTION  FUNCTION 


The  integral  of  Eq.(49)  is  a  function  found  in  other  applications  than  the  one 
discussed  in  this  paper.  It  is  desirabla  to  have  graphs  of  this  function  available 
in  general  form  rather  than  the  specialised  form  of  Fign. 1-50,  No. 1.  The  integral 
is  s  special  csss  of  the  incomplete  Toronto  function*  described  by  Heatley<4a>  and 
Fishar<1T),  which  is  defined  as 


r) 


tm’ner‘aIn(2rt)dt  . 


tt 

(100b) 


Using  this  notation,  Eq.(49)  for  the  cumulative  distribution  function  becomes 


Pj,  •  1  -  Tyy^  { 2N- 1 , N- 1 , /Nx) 


(100c) 


*  la  normal  correlation  tlioory,  ths  quantity 


?(V 


ia  (iron  by  Fishor'**)  aa  the  Halting  lore,  (or  largo  aaaplas,  of  tha  frequency  element 
of  tha  quantity  B *  ■  a  whara  A*  it  tha  staple  aatieata  of  tha  auitipla  correlation 
cooffieiant  of  a  random  variable  y  with  other  variabiea  *l,  «,••••  s  ,  n,  ia  the  eiie 
of  tha  eaaple,  end  4g  *  n» fia  where  p  ia  tha  population  Multiple  correlation  coefficient. 
The  cueuletive  dietnbution  ie 


end  can  ba  obtained  froa  the  curvea  in  Figa<13  to  32. 


The  function  plotted  in  Figs.  13  to  32  is 
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The  distribution  function  for  the  sum  of  two  variates  of  noise  alone  is  ob¬ 
tainable  by  use  of  the  convolution  theorem,  giving 


:pj  "  1  *78'  MVl)er/  7 


(105) 


and  the  cumulative  distribution  is  also  obtainable,  giving 


f(y)<*y 


’4  ^  - 4  ,v 

me  *  T~  ye  4  erf  -*• 


(106) 


However,  when  (V  >  2  there  seems  to  be  no  closed  solution  corresponding  to  Eq.(lOS) 
or  (106),  Since  these  cases  are  for  noise  alone,  the  signal-plus-noise  situation 
must  be  attacked  by  other  means. 

It  turns  out  that  if  the  momenta  of  the  distribution  for  one  variate  are  known, 
the  moments  of  the  distribution  for  the  sun  of  N  variates  may  be  found  directly. 
Formulae  are  given,  for  instance,  in  Cramer,  page  345, <10)  for  the  first  few  central 
momenta,  which  are 


Mj  ■ 

^4  ■  (Vm4  ♦  ZNiN-Dp* 

Mg  %  10 (VVJ  ♦  15/VVJ  *  15(VV2M4  -  45iVVJ 


The  corresponding  coefficients  in  the  Gram-Charlier  series  then  become 


C3  “  -3JiV/a 


(107a) 


(107b) 


(107c) 


( 107d) 


Cfl  "  6  \N  ’ 


(108a) 


(108b) 


(108c) 
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The  a1 s  in  Eqs.(108a-c)  are  the  central  standard  moments  for  one  variate.  Note 
that,  in  the  square  law  case,  if  N  is  put  equal  to  1  in  Eqs.(85c-f)  and  the  resulting 
M's  used  in  Eqs. (107a-d),  the  /x' a  for  N  variates  are  correctly  given.  If  a  moment 
generating  function  can  be  found  for  the  case  of  N  variates,  then  it  is  immaterial 
which  method  is  used;  but  in  the  case  in  which  such  a  function  is  not  available,  the 
Eqs. (107)  must  be  used  (or  some  method  essentially  equivalent). 

To  handle  the  linear  detector  it  is  now  sufficient  to  find  the  moments  for  one 
variate  only.  Rice*1**,  page  107,  gives  the  required  expression  as 


"  *i/3r(ui)iFi  (-2*1'-*)  • 


(109) 


Rice  also  gives  the  first  two  moments  as 


(110a) 


va  ■  2(1+*) 


(UOb) 


$ 

i’fi 

W 

«! 


To  calculate  one  needa  to  know  the  function, ^(-3/2,1,-*).  Thia  may  be  obtained 
by  use  of  the  recursion  relation 


ijf, (<*♦!,  c,  z)  ♦  (o-c),Fj(o-l  ,c,i)  ■  (2n  +  z-c)  ,F,(a,c,  z)  (111) 


•  w 

v/ 

% 

U 


by  putting  a  «  ,  c  ■  l,  z  «  -x.  The  result  is 


2J/j  (2  +  *)  -Vf*"7lo(f) 


(110c) 


V. 


■  i(2*ix*x2)  . 


( 1 10  d ) 


II 


**  \«  »>•»>  ^  yw>  >>  v>  w  uvvxyxuvwuN  (.vwi/K'Wvvy  '.wkww.m 


This  corresponding  central  moments  are 


Ma  ■  <r3  ■  2(l+x)  -  v3 


-  2v3  -  2v,(l  +  2*)  -^ye‘7I0(y) 

/x4  -  4(2+4*+*J)  -  3v*  -  4v,  [d-*)v,  4(f)] 


(112a) 


(112b) 


(112c) 


The  standard  central  moments,  and  then  the  c’s  of  Eqs. (108a-c) ,  are  directly  obtain¬ 
able  from  these  formulae,  though  the  process  is  somewhat  tedious  due  to  the  cumbersome 
form  of  Eqs. (112a-c) .  The  functions  v,  to  v4  are  shown  graphically  as  a  function  of 
x  in  Fig. 34. 

To  obtain  the  bias  level  Yb  for  the  linear  detector  for  N  >  2,  one  can  use  the 
G.C.  series  for  noise  alone.  Setting  *  *  0  and  v,  ■  V'n/i  in  Eqs.  (,112a-c)  gives 


and 


At,  2  -  -  0. 429  (113a) 

/u.j  “V?( 7T”3)  “  9*1?72  (113b) 

£i4  *  8  -  ■  0.598  (  113c) 
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The  cumulative 


distribution  function  is  now  equated  to  IjJ ,  giving 


0.  693 JV 
n 


0. 1053 
*  JVVJ 


d>2(T) 


0.0108 
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0.00555 

JV 


0*(T) 


--•(116.) 


where 


T 


For  any  given  n  and  N,  Ffcmay  be  found  from  Eq.(116)  by  trial  and  error  methods. 
If  an  approximate  value  of  T  is  found  by  neglecting  all  but  the  first  term  in  Eq.(116), 
a  more  accurate  value  obtained  by  Newton’s  method  is 


(117) 


It  is  better,  however,  to  plot  Eq.(116)  giving  n  as  a  function  of  T  and  N  from 
which  is  finally  obtained  the  bias  level  graph  of  Fig.  35  showing  Yb  as  a  function 
of  n  and  N  for  the  linear  detector. 

Since  for  finding  the  bias  level  it  is  necessary  to  know  the  distribution 
functions  only  for  large  values  of  the  argument,  it  is  possible  to  find  an  ap¬ 
proximate  solution  valid  in  this  region.  Consider  a  distribution  function  given  by 


dP  ■  ve  "3"  dv 


(117a) 


for  v  going  from  -oototoo.  The  Nth  convolution  of  this  function  will  be  nearly 
the  same  for  large  values  as  if  (117a)  went  only  from  0  to  ®,  because  the  large 
values  in  the  sum  of  N  variates  are  most  probably  produced  by  addition  of  large 
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values  of  every  variate,  and  for  large  values  (in  fa  fo’*  :.e  values) 
the  two  distribution  functions  are  identical.  The  characteristic  function  of  Eq. 
(117a)  is  given  by  pair  710.1  of  Campbell  and  Foster  to  be 


C 


(U7b) 


For  the  sum  of  N  variates 


£  Xtl 

CN  •  (-1)*(2 rr)Ve  S 


(117c) 


The  probability  density  function  is  obtained  from  pair  740.2  of  Campbell  and  Foster  as 


y  >>  1 


( 117d) 


where  D #  is  the  parabolic  cylinder  function  of  order  N.  In  terms  of  the  derivative 
of  the  error  integral  as  defined  in  Eq.(62), 


y>>  1 


( 117e) 


Note  that  for  N  ■  2,  Eq.(117e)  becomes 


( 1 17  f ) 


Referring  to  Eq.(lOS),  the  exact  expression  for  this  case,  it  is  seen  that  Eq.(117f) 
can  be  obtained  by  neglecting  the  first  term  and  replacing  srf  y/2  by  1,  both  of 
these  approximations  being  very  good  if  y  »  1. 


The  approximate  cumulative  distribution  is  easily  obtained  from  Eq.(117e)  by 
direct  integration  and  gives 


( 117g) 
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The  bias  level  is  easily  obtainable  from  this  expression  by  equating  it  to  Q  and 
solving  for  ^  by  means  of  the  tables  of  <t> *,  or  by  plotting  graphs.  The  method  is 
not  very  practical  for  N  >  20  since  suitable  tables  do  not  exist. 

It  is  interesting  to  note  that  no  such  approximation  as  Eq. (117g)  is  obtainable 
for  the  square  law  case. 

Graphs  of  the  probability  density  functions  for  signal -plus-noise  have  been  ob¬ 
tained  by  numerical  convolution  for  some  selected  cases  and  are  shown  in  Figs, 36  to  41. 


RESULTS  OF  THE  LINEAR  DETECTOR  CALCULATIONS 

The  difference  in  results  for  the  linear  and  square  law  detectors  turns  out 
to  be  so  small  that  extreme  accuracy  must  be  used  in  the  calculations  to  show  the 
relation  in  its  true  form.  One  such  comparison  graph  was  calculated  and  is  shown 
in  Fig. 42.  Also,  in  Fig.  43  is  shown  the  difference  in  db  in  the  two  cases  at  P  • 
0.50.  The  two  are  identical  at  N  •  1,  the  linear  law  becomes  better  by  a  maximum 
of  0.11  db  at  N  ■  10,  the  two  are  again  equal  at  N  •  70,  and  the  square  law  then 
becomes  better  and  asymptotically  exceeds  the  linear  law  by  0.19  db  as  N  ®  having 
reached  0.16  db  at  N  ■  10C0.  These  results  show  conclusively  that  there  is  little 
to  choose  between  the  linear  law  and  square  law  as  far  as  theoretical  signal  threshold 
is  concerned. 


EXPANSIONS  IN  LAGUERRE  SERIES 

In  certain  cases,  particularly  for  low  values  of  N,  the  Gram-Charlier  series 
may  not  be  the  best-suited  type  of  expansion  for  distribution  functions  which  are 
zero  for  all  negative  values  of  the  amplitude.  For  low  values  of  N,  a  suitable 
expansion  for  such  functions  is  the  following: 


(118) 

1*0 

where  L“(y)  is  the  generalized  Laguerre  polynomial  defined  by  Eq.(82),  or  by 


LJ 


e*  z'4  d' 


3? <*■**"■> 


(119) 


The  orthogonality  relation  which  makes  the  expansion  possible  is 


± 1_  r 

a+i*  1  )  J 


r(a+t  +  l )  I  ®‘*  za^(z)L;a(z)  -  Si;. 


(120) 
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Since  there  are  two  arbitrary  constants,  a  and  (3,  in  the  expansion  of  Eq.(122),  it 
is  possible  to  make  c{  ■  ■  0  by  a  proper  choice  of  a  and  j3.  These  relations  are 
easily  determined  by  equating  Eqs.(125)  and  (126)  to  zero  and  solving  simultaneously. 
The  results  are 


and 


(127) 


(128) 


1 

G°  '  /8r(a+l) 


(129) 


(130) 


1 

/3r(a*4) 


(a+3 ) 


(131) 


The  coefficients  past  cQ  are  so  complicated  that  the  whole  value  of  this  type  of 
series  seems  to  depend  on  the  fact  that  the  first  term  alone  is  often  a  good  ap¬ 
proximation.  This  approximation  is 


dP 


(132) 


and  the  corresponding  cumulative  distribution  function  is 


P 


/%  m 

v} 

1  f(y)  dy  •  1  -  I 

l.-L  -i 

K 

f  a1 

(133) 
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where  I  is  the  incomplete  gairma  function  as  defined  by  Eq.(40). 


There  is  a  striking  analogy  between  Eq.(132)  and  the  corresponding  normal 
approximation.  In  both  cases  the  distribution  for  the  sum  of  N  variates  is  simply 
obtained  by  multiplying  both  and  <r3  by  N,  A a  iV  -*  ®,  both  the  normal  approximation 
and  Eq.(132)  approach  the  true  distribution  (and  each  other).  In  any  particular 
case,  however,  the  convergence  properties  of  one  approximation  will  be  more  useful 
than  the  other. 

In  the  square  law  case,  for  *  ■  0,  vt  ■  N  and  cr3  ■  N.  Substitution  of  these 
values  in  Eq.(132)  gives 


dP  ■  7TTT  e"*y*"‘dy  .  (134) 


Note  that  this  is  the  same  as  Eq. (38),  the  exact  expression.  Thus  in  this  particular 
case  the  first  term  gives  the  whole  correct  result.  The  third  coefficient  from 
Eq.(131)  is  easily  shown  to  be  zero,  as  all  the  following  coefficients  will  be. 

j  In  the  square  law  case  where  *  ^  0,  v,  ■  N(  1+*)  and  era  ■  N{ 1+2*).  Substitution 
of  these  values  in  Eq.(132)  gives 


dP  - 


(135) 


and  from  Eq. (133) , 
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1 


I 


----- * - 

/\( l>2x) 


/V(  1+s)3 
1  +  2* 


1 


(136) 


A  comparison  of  the  particular  case  N  ■  3,  x  ■  1  is  shown  in  Fig. 44.  Curves  are 
given  for  the  exact  distribution  function  (Eql37))  and  the  two  approximations  given 
by  Eqs. (96)  and  (135). 

For  the  linear  case  with  x  ■  0,  v,  ■  /V/rr/2  and  <r3  ■  N(2-^/2),  the  cumulative 
distribution  is,  from  Eq.(133), 
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1  ~  I 


//V(2-tt/2  )  '  4/tt-i 


-1 


(137) 
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OTHER  SERIES  APPROXIMATIONS 
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It  is  theoretically  possible  to  develop  still  other  series  approximations  for 
the  various  distribution  functions.  Forainatancel  it  might  be  thought  advantageous 
to  use  a  sum  of  terms  of  the  type  yae**  particularly  in  the  linear  case.  While 
this  turns  out  to  be  possible,  even  the  first  coefficient  is  so  difficult  to  calcu¬ 
late  that  the  process  is  impractical. 


METHODS  OF  INTEGRATION  INVOLVING 
SUBTRACTION  OF  NOISE 

Certain  practical  difficulties  arise  in  maintaining  the  bias  level  at  the 
correct  value  in  an  electronic  detector,  particularly  if  the  number  of  pulses  inte¬ 
grated  is  large.  The  trouble  may  arise  from  fluctuations  in  amplifier  gain,  the 
bias  supply,  or  the  noise  level  itself. 

A  solution  of  this  problem  is  to  have  the  gain  of  the  amplifier,  or  the  bias 
level,  or  both,  controlled  by  some  sort  of  average  value  of  the  noise  output.  Ob¬ 
viously  the  time  constant  of  the  control  device  must  be  neither  too  long  nor  too 
short.  One  scheme  which  has  been  used  is  to  subtract  a  pulse  known  to  consist  of 
noise  only  from  each  possible  aignal-plua-noise  pulse*  (see  paragraph  3,  paga  11, 
No.l).  Thus  in  the  absence  of  a  signal ,  the  average  value  of  any  number  of  composite 
pulses  will  always  be  sero,  and  the  required  bias  level  will  be  comparatively  low. 


DISTRIBUTION  FUNCTIONS  FOR  COMPOSITE  PULSES 
OF  SIGNAL-PLUS-NOISE  MINUS  NOISE 

When  a  noise  pulse  is  subtracted  from  each  signal -plus-noise  pulse,  the  theo¬ 
retical  distribution  functions  will  be  entirely  different  from  previous  cases.  The 
square  law  case  is  the  only  one  that  can  be  treated  in  any  reasonable  fashion.  The 
distribution  function  for  one  variate  of  signal  plus  noise  is  given  by 

dP  -  /0(2V7P)  dY  (138) 


and  the  characteristic  function  is 


C 


(139) 


*  This  subtraction  con  b*  accomplished  by  nssna  of  s  gste  which  operates  st  double  ths 
repetition  frequency.  0»  every  other  gate  only  a  noise  pulae  of  reversed  phase  goes  through 
the  integrator. 
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Subtracting  a  positive  noise  variate  is  equivalent  to  adding  a  negative  noise  variate. 
The  distribution  function  for  a  negative  noise  variate  is 

dP  -  b*  Y  <  0  (140) 

■  0  Y  >  0 


and 


C 


(141) 


To  obtain  the  characteristic  function  for  the  sum  of  a  variate  from  the  distributions 
of  Eqs.(138)  and  (140)  it  is  only  necessary  to  take  the  product  of  the  characteristic 
functions  given  by  Eqs.(139)  and  ( 141) »  giving 


C  - 


(142) 


This  is  the  characteristic  function  for  one  so-called  coeposite  pulse.  The  characteris¬ 
tic  function  for  the  sum  of  N  composite  pulses  is  simply 
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(1 -p2)n 


(143) 


In  the  case  of  noiae  alone  (x  «  0), 


C 


1 

(1 -p*)" 


(144) 


and  the  anticharacteristic  function  is,  by  pair  569  Campbell  and  Foster, 


dps 


1 

v^(/v- 1)! 


Y  |  dY 


(145) 


where  Kv _ ^  is  a  modified  Bessel  function  of  the  second  kind  and  is  given  by  the 
finite  series 
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The  cumulative  distribution  for  the  sum  of  N  composite  noise  variates  may  be  found 
by  use  of  the  series  (146)  and  term  by  term  integration.  However,  for  (V  greater 
than  3  or  4  che  process  rapidly  becomes  impractical. 

Again,  it  is  necessary  to  find  momenta  and  proceed  by  meana  of  Gram-Charlier 
series.  For  noise  alone,  the  moments  are  easily  found  from  Eq.(l44)  to  be 


(147) 

,  i  ever 
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in  particular, 


■  2/V  ■  cra 
M4  •  12N(N*1) 


a 


« 


15  ♦ 


45  30 
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The  only  coefficients  different  from  zero  in  the  first,  six  are  c  and  c  to  the 
order  of  1/iV, 


Thus 


(154) 


(155) 
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and  similar  Co  Eq. (98)  is  the  cumulative  distribution 


(156) 


The  bias  number  is  found  by  setting  this  expression  equal  to  and  plotting  Yk 
as  a  function  of  n  and  N.  Results  are  given  in  Fig.  45.  In  the  special  ease  N  ■  1, 
the  cumulative  distribution  function  ia  simply  e'r/i  for  Y  >  0,  and  the  bias  number 
is  obtained  from  this  expression  rather  than  from  Eq.(156).  The  anticharacteristic 
function  of  the  general  case,  Eq.(143),  may  be  obtained  by  use  of  the  convolution 
theoremt  pair  202 »  Campbell  and  Foster.  Let 
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Then  from  Eq. (37), 
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and  by  pair  525. 2,  Campbell  and  Foster, 
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(160) 


Applying  the  convolution  theorem  gives 

N-\ 

-Y.Nm  /**  2  ___ 

df"'  ‘"’(Fin  J, (/£)  i-’o  (i6D 

For  Y  <  0,  the  lower  limit  of  the  integral  in  Eq.(161)  is  0  rather  than  Y, 
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PROBABILITY  DENSITY  FUNCTIONS  FOR  Y  <  0 


To  evaluate  the  integral  in  Eq.(161)  when  the  lower  limit  is  zero  is  straight¬ 
forward  but  tedious.  First  one  evaluates  the  integral 


/(*> 


e-s*  4.,  (2 /Nxy)dy 


(162) 


by  use  of  characteristic  functions  in  a  manner  entirely  similar  to  that  used  in 
Eqs.(74)  to  (80).  The  characteristic  function  of  the  function  of  Eq.(162),  with 
*  *  0,  is 
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and 
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Then  by  expanding  one  obtains  the  coefficient  of 
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(N-i)  ?  (-n*-1-* 

k]  (N-l-k)l 


(165) 


and  from  Eqs. (164),  (165),  and  (161), 
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In  terms  of  Laguerre  polynomials,  using  Eq. (82) , 
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The  first  few  polynomials  are  given  in Eqs. (124a-d>.  Some  special  cases  of  Eq.(l68)  are 


e* 

N  -  1  dPt  •  Y  <0  (169a) 

e *  /  \ 

iV  -2  dpj  -  «iy— (i+f-yj  y<o  (169b) 


N  »  3 


dP, 


y<o  . 


(169c) 


The  cumulative  distributions  for  Y  <  0  may  easily  be  obtained  by  integrating  (l69a-c). 
Obviously,  the  expressions  in  Eqs. (167) and  (168)  are  practically  useful  only  for 
small  values  of  (V. 


PROBABILITY  DENSITY 
FUNCTIONS  FOR  Y  >  0 

To  find  a  general  expression  for  Eq.(161)  giving  the  distribution  function 
when  y  >  0  is  a  task  of  tremendous  proportions.  Consider,  for  instance,  the  special 
case  /V  ■  1,  Equation  (161)  becomes 


dPj  -  dYe 


e‘3r  I0(2/xy)dy 


Y  >  0  . 


(170) 


By  means  of  the  substitution  y  ■  vJ/4,  this  becomes 


K"f  /v®  v3** 

dPj  -  dYS-^-l  e'^I^v/Tjdv  Y >  0  .  (171) 
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This  can  be  expressed  in  terms  of  the  Q  function  defined  by  Eq.(16). 


Eq.(169a)  was 


dP, 


dY*-r-Q{/x,  2/7)  Y>  0 


dPt  •  Zj-dY-  Y<  0  .  , 


(  172) 
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Thus,  for  the  case  /V  ■  1,  the  whole  distribution  function  is  described  by  Eqs.(169a) 
and  (172).  A  graph  of  this  function  for  various  values  of  *  is  shown  in  Fig. 46. 
Note  that  if  x  *•  0,  Q(0,2/F)  «  e"2*',  and  Eq.(172)  for  Y  >  0  reduces  to 


e-y 

dP,  •  -j-dY  }'>0 

and  from  Eq. ( 169a) 
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when  *  •  0.  Thus  over  the  whole  range  of  Y 


(175) 


which  checks  Eq.(145)  when  N  ■  1. 
For  N  ■  2,  Eq.(161)  becomes 
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dYcy-uJ  (^)7(y-F)e*Jy  I,(2/J7y)dy 
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This  integral  may  also  be  evaluated  in  terms  of  the  Q  function.  The  process  requires 
a  large  number  of  integrations  by  parts  and  is  very  time-consuming.  The  result 
turns  out  to  be 


dP, 


(i+Z-y'jQiSzZ,  2/7) 


This  equation  is  already  so  complicated  as  to  be  nearly  useless.  Thus  it  was  not 
thought  worth  while  to  seek  a  general  expression  of  this  type  for  arbitrary  N  when 
Y  >  0. 

Note:  If  r  *  0  in  Eq.(77),  it  reduces  to 

S~Y 

dpt  *  dr-j-u+n  r>o  (178) 


which  may  also  be  obtained  from  Eq.(169b)  by  substituting  -T  for  Y. 


CUMULATIVE  DISTRIBUTION  FUNCTIONS 

The  effort  in  Eqs.(161)  to  (178)  has  been  concerned  with  obtaining  the  proba¬ 
bility  density  function  for  N  variates  of  signal-plus-noise  minus  noise.  To  find 
the  cumulative  distribution  functions  exactly  is  difficult,  especially  for  Y positive. 

A  case  which  can  be  solved,  however,  is  that  for  N  •  1.  For  Y  negative  the 
answer  is  simply  obtained  from  Eq.(169a)  and  is 


/>  -  1  -  Y  <0  .  (178a) 


For  Y  positive,  using  the  result  of  Eq.(172), 
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e*Q(/T,2/7)dy  Y>  0  . 
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Since  the  value  of  at  Y  ■  0  is,  from  Eq.(178a),  1  -  e  */2,  Eq.(178b)  may  be  re¬ 
written  as 
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but  from  the  definition  of  Q  in  Eq.(16), 
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Replacing  Q  by  its  defining  integral  in  Eq.(178c)  gives 


P,  ■  1  -  ^  e“s,I0(  2*^x7  )dz  . 

(178e) 

Integration  by  parts  is  now  used,  letting 

u  ■  J*  e~2*  I6(2/xi)dz  | 

dv  ■  dy  j 

►  ( 178  f) 

da  ■  -e‘2yI0(2/jty)dy  | 

v  ■  e*  t 

/*  * 

\uv\*  •  er  J  e~2lI0(2/xy)ds  -  ~ 

(178g) 

7 

-  -  T  • 

(178h) 

Thus 

u  X 

e  1  rY 

P,  ■  1  -  —  Q(/T,2/y)  +  e~*  1  vdu 

2  Jo 

(178i) 

or 

/“I  /•* 

P,  -  1  -  —  Q(/T,2/?)  -  e*  J  e->I0(2Sxy)dy  . 

Jo 

( 17 8 j ) 

The  integral  term  in  Eq,(178j)  is  just  1  -  Q(/ Tx,/2?\  and  the  final  result  is 

I'-f 

P,  -  <?(/2l,/2?)  -  —<?(/*,  2/?)  y>0  . 

(178k) 

a 

For  x  ■  0,  Q(0,j3)  •  e  1  ,  ard 

P,  -  £ 

(1781) 
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agreeing,  as  it  should,  with  the  result  obtained  from  Eq.(145)  by  letting  N  ■  1 
and  integrating.  For  Y  ■  0,  Q(a,0)  ■  1,  and 


Px 


(178m) 


agreeing  with  Eq.(178a)  when  Y  ■  0. 

The  bias  number  for  use  with  Eq.  (1781c)  is  obtained  by 


r* 


0.693 

n 


2 


or 


■  2.30  log,0  n  -  0. 327  . 


(178n) 


(178o) 


In  Fig.  47  is  shown  a  graph  comparing  Eq. (178k)  with  Eq.(23)  for  n  ■  10*,  where 
P  ia  plotted  as  a  function  of  *. 

Though  it  might  be  possible  to  calculate  the  cumulative  distributions  for  N 
>  1  by  a  method  similar  to  that  uaed  for  N  •  1,  it  would  be  very  tedious.  Therefore 
resort  ia  made  to  Gram-Charlier  aeries,  as  before.  The  momenta  are  directly  obtainable 
from  the  characteristic  function  given  in  Eq.(143), 


JL 

V, 

dp' 

p  ■  o 


(179) 


There  seems  to  be  no  readily  obtainable  expression  for  1/  in  closed  form.  The  first 
six  moments  obtained  directly  from  Eq.(179)  are;  * 


vi  "  (180a) 

■  ( Nx)a  +  2 Nx  +  2JY  (180b) 

-  (Nx)3  +  6 (Nx)a  +  6Nx(N*l)  (180c) 

v4  -  (Nx)*  *  12(Nx)3  *  12(iVx)3(iV+3 )  ♦  24(1¥+1)  +  \2N(N*l)  (180d) 


v,  -  (Nx)s  *  20 (Nx)*  +  20((Vx)3(iV+6)  +  120  (Nx)3  (iV+2 )  ♦  60iV*(iV+l )  (JV+2 )  . 

( 180e) 
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Vt  -  OV*)4  ♦  30  (iV*)5  ♦  3Q(iV*)4  (JV+10)  ♦  1200Y*)3(3(V+10) 
+  180OY*)3  (/V+l)  (iV+6)  +  360JV*'JV+l)((V+2) 

+  120/V(iV+l)(iV+2)  . 


( 180  f ) 


The  corresponding  central  moments  are: 

a 

■  2JV*  ♦  2JV  "  2iV(  l  ♦* )  «  a 
Mj  ■  6jV* 

M4  -  12(iV*)J  ♦  24/V*(/V+l)  ♦  12N(N*1) 


(181a) 

(181b) 

(181c) 


M,  -  120((V*)3  ♦  360 (/V*)J(/V+3)  ♦  360(V*(:V+1)  (:V+2)  ♦  120A((Y*l)(/V+2)  • 

(laid) 

The  central  standard  moments  are: 


a,  - 


3* 


/2N(1**); 


(182a) 


3  ♦ 


3 ( 1*2*) 

N{  l+*)J 


(182b) 


45 (3x2*ix*l) 


?y  15  ^  ...  1 1 

4  A(  1  +  *)3 


(182c) 


The  coefficients  of  the  series  are,  from  Eq.(71): 

x 

C3  "  "  1 

2/2A(l+*)3 

1*2* 

C<  “  8A(1  +  *)3 

s3 

C*  "  16(V(  1 1*)3 


(183a) 


(183b) 


( 183c) 
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The  Gram-Charlier  aeries  for  the  probability  density  function  is  given  by  Eq.(93)  where 


y-^i  . 

t  ■  — ■£— ,  t'j  ■  Sx,  a  ■  /2/V(  1  +  Jt) 


(184) 


and  the  cumulative  distribution  is  given  by  Eqs.(98)  and  (99).  Figures  53  and  54, 
No.l,  showing  the  comparison  between  the  ordinary  case  and  the  composite  case,  were 
computed  using  the  Gram-Charlier  series  developed  above.  There  appears  to  be  no 
significant  difference  in  the  probabilities  of  detection  for  N  between  1  and  10.  For 
N  between  100  and  1,000,  the  composite  case  gives  an  effective  signal-to-noise  ratio 
about  1  db  lower  than  the  ordinary  case. 


ANOTHER  APPROACH  TO  THE  DETECTION 
CRITERIA- PROBABILITY  THAT  SIGNAL- 
PLUS-NOISE  EXCEEDS  NOISE  ALONE 

The  method  of  setting  a  bias  level  and  calling  any  aignal-plus-noise  or  noise 
alone  which  exceeds  this  level  a  signal  is  not  the  only  possible  way  of  defining 
detection.  Another  method  is  based  on  asking  what  is  the  probability  that  any  given 
signal  will  be  larger  than  any  noise  pulse  during  a  given  interval  of  time'1**.  The 
interval  of  time  taken  would  logically  be  the  false  alarm  time,  as  defined  previously. 
In  this  time  there  will  be  n/N  ■  n'  independent  groups  of  noise  pulses.  If  the  proba¬ 
bility  that  a  single  integrated  group  of  signal-plus-noise  pulses  exceeds  a  single 
group  of  noise  pulses  is  called  /J(*,A(),  then  the  probability  that  the  group  of 
signal-plus-noise  pulses  exceeds  all  of  the  n'  groups  of  noise  pulses  is  simply 

P  -  (P.U.N)]"'  (185) 


This  probability  is  a  little  difficult  to  interpret  properly.  It  means  that  if 
during  the  false  alarm  time  a  signal  of  strength  x  appears,  it  will  have  this  proba¬ 
bility  of  being  larger  than  any  noise  pulse  group  appearing  during  the  same  time.  The 
difficulty  is  how  to  pick  out  the  largest  signal  over  a  period  of  time,  and  what 
to  do  when  many  signals  are  present.  These  are  reasons  why  the  earlier  detection 
criteria  are  thought  to  be  superior,  since  they  provide  clear  answers  for  the  above 
questions.  The  criteria  presented  above  may  be  of  special  value,  however,  when  a 
target  is  known  to  be  present.  Such  is  the  case  when  a  target  is  being  automatically 
tracked,  and  one  wishes  to  calculate  the  probability  that  it  will  be  subsequently 
lost  due  to  the  noise  exceeding  the  signal. 

The  probability  density  function  for  A'  signal-plus-noise  pulses  minus  N  noise 
pulses  has  been  indicated  in  Eq .(161). 

To  obtain  the  probability  that  the  sum  of  iV  signal-plus-noise  pulses  will  be 
greater  than  S  noise  pulses  it  is  only  necessary  to  integrate  Eq.(161)  from  0  to 
°°.  It  will  be  easier  to  obtain  the  probability  that  /V  noise  pulses  exceed  N  signal- 
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plus-noise  pulses,  however,  since  this  requires  the  integral  from  -®  to  0,  and  an 
expression  is  available  for  Y  <  0  in  Eq.(167).  Thus 


/•<»  kns-i  (,V*A-1)! .F.(-k,  N, 

.  -fJ-  I  y  - ‘7,  2/(-na-‘  (186) 

(N-l) !  U  (iS-k-l)\k\2N*h 


Now  one  substitutes  z  for  -Y  and  interchanges  the  summation  and  integration  signs, 
obtaining 


p  ‘  (iv-D!  2-  (iV-*-l)!*!2A'+‘  Jo 


Hie  integral  is  simply  (JV-fc- 1)!,  and  therefore 


_  e  r  (iV»*-l)!  Ff  ,  „  1 

2  S  (N-l)!*!2^-“  •(  '  *“ 


Or  in  terms  of  Laguerre  polynomials,  using  Eq.(82), 


it* 

'T  *■#-» 


From  Eq.(188)  a  more  convenient  form  may  be  obtained  by  introducing  a  dumny  index 
i  and  interchanging  summation  signs,  leading  eventually  to 


N*  fNx\ * 

.  £Z  V*  W  (iV^-1)? 

2  2-  if  (iVwi-1)!  2f  (*-i)!2Ar+fc’‘ 


The  outside  summation  in  Eq,(190)  is  obviously  a  polynomial  in  x  of  the  N-lth  degree 
and  with  N  terms.  It  is  rather  curious  to  note  that  if  one  puts  i  ■  0  in  Eq.(190), 
the  following  identity  results: 


2*-t  ,  *^1  (MiliL 

2w  (N-l)  !*!  2J 


VjBM  V  M±\  HWI.SUJ  Ml.-  I  SL/I 


In  ocher  words,  the  constant  term  in  the  polynomial  is  always  unity. 


The  first  few  cases  for  low  values  of  N  are: 


p,  ■  A  «  * 


4 


5* 


(192a) 

j) 

(192b) 

9s  9xJ\ 

16*128  / 

(192c) 

29s  x3  Xs  \ 

*  32  *  4  *48  ) 

(I92d) 

325s  57  5s*  4375s*  625s4  \ 

*  256  *  1024  *  6144  '  32,  768  /  * 

(192e) 

Obviously  for  N  very  large,  these  expressions  rapidly  become  useless,  and  it  is 
necessary  to  use  the  Gram-Charlier  series  of  Eqs.(184)  and  (98).  The  lower  limit 
Y  is  replaced  by  zero,  giving  for  the  series 


P  *  i[w‘m]  ♦  c3<t>HT)  -  c^3(T)  -  ct4>s(  T)  -■ 


(193) 


where 


v»  ,/T" 
T  ~ 


(194) 


and  Cj ,  c4  and  c4  are  given  by  Eqs. (183a-c) .  A  graph  of  P  as  a  function  of  s  and  N 
is  shown  in  Fig. 48.  For  very  small  values  of  P,  more  terms  may  be  necessary  in  the 
series  of  Eq. (193). 


USE  OF  CUMULANTS  IN  OBTAINING 
GRAM-CHARLIER  SERIES  COEFFICIENTS 

It  is  often  much  simpler  to  obtain  the  cumulants  for  a  given  distribution 
function  rather  than  the  various  moments.  The  cumulants  may  be  defined  by 


l08<  c),.. 


;  195) 


-*r_  •"  L "21 ^ 7 S' .  **J" -T-**-'  n  ■» , x  \  ^J* a«'v-.  v «s.T rsS7 rwT sLTmJT mD T iwT m  mm'm'm  m-f. 
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where  C  is  the  characteristic  function  of  the  given  probability  density  function 
(see  pages  61-65  of  Kendall<#)).  The  cumulants,  except  the  first,  are  invaria-t 
with  respect  to  a  change  of  origin.  Also,  for  the  distribution  of  the  sum  oi  N 
variates,  it  is  only  necessary  to  multiply  every  cumulant  by  N,  as  is  evident  from 
the  defining  Eq . ( 195 ) .  The  coefficients  of  the  Gram-Charlier  series  in  terms  of 
cumulants  are  given  on  page  149  of  Kendall.  The  cumulants  in  standard  measure  may 
be  defined  as 


K, 


(196) 


In  terms  of  standard  cumulants,  the  coefficients  of  the  series  are: 


c0  “  1  cl  '  C2  “  0 


K, 

C»  “  3! 
“  4! 


(197a) 

(197b) 

(197c) 

( 197  d ) 


The  first  term  in  Eq.(l97d),  K#,  is  omitted  in  the  0, 3,4,6  approximation. 
Consider  the  square  law  case  where,  from  Eq.(35), 

«  * 

c  e'*e^ 

C  *  "  . - 

P  +  1 


From  Eq. ( 195) , 


For  N  variates, 


and 


log.  C  »  -  x  +  -  -  In  (p+l) 

p  +  l 


^■(i-OJtis  +  l)  if  1 


k.  -  N(  i- 1 ) !  (  ix+ 1 )  . 


*i  (  i-l)  !(  t*+l) 
K.  ■  — •  ■  • 


(198) 

(199) 

(200) 

(201) 

(202) 
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In  particular, 


2( 3*+l) 


NUi(  2 


*+l)^ 


K 


6<4x+l) 
4  "  JV(2*  +  1)J 


(203a) 


(203b) 


and  it  is  at  once  evident  that  cJ(  c4,  and  cf  obtained  from  Eqs.(197b-d)  are  iden- 
tical  to  the  valuea  given  by  Eqa.(92c-e)  by  means  of  a  much  longer  process. 

In  the  case  of  a  composite  pulse  of  signal*plus»noise  minus  noise,  the  charac* 
teristic  function  is  given  by  Eq.(142)  and 


log  C  ■  -  *  +  — —  -  log  (l -p3) 

*  p*  1  * 


(204) 


Again  by  means  of  Eq.(195)  it  is  easy  to  derive,  for  N  variates, 


Ki  ■  iV(  i-1) !  [ix+2]  i  even 

•  N(  i-1) !  ( ix)  i  odd,  /  l 


(20S) 


or 


and 


-  iV(  t-l) !  [ix+l  +  t-l)*] 


(206) 


Special  cases  are: 


„  *i  (i-0![i*  +  lM-l)‘l 


(207) 


fix 


Nl/1  [2  ( *  + 1  )J 


(208a) 


3(  2xf  1 ) 
Nix*  l)2 


(208b) 


and  again  by  Eqs. ( 197b-d)  the  coefficients  are  seen  to  be  the  same  as  given  by 
Eqs. ( 183a-c) . 

S4 


L 


\WWV\  vA  V\ 

w  w’m  awsweusti  sea*  m-u  ■  v.  arw  iru  Srw  m’umv;  rfvi  a*v  e#Vj  arvj  arv<  arvj  rfijjrv  ai-u  wvuni  1 


In  a  cftfte  such  as  the  linear  one  where  the  characteristic  function  cannot  be 
obtained,  the  cumulants  are  still  useful  and  may  be  found  from  the  moments  by 
means  of  the  formulae  at  the  bottom  of  page  63  of  Kendall.  The  first  few  are; 


(209a) 


(209b) 


K 


i 


Vi  "  3VJV1  +  2Vl 


(209c) 


are 


■  v4  -  *VjV,  -  +  ISVjVj*  -  6v*  .  (209d) 

Tl»  are  now  obtained  by  multiplying  by  N  and  dividing  by  o-‘.  The  coefficients 
then  obtained  as  before  by  Eqs. (197e»d). 


BEST  POSSIBLE  DETECTOR  LAW 

It  is  of  considerable  importance  to  know  whether  there  may  be  some  detector 
law  which  will  give  results  which  are  appreciably  better  than  the  linear  or  square 
law  cases  which  have  already  been  considered. 

The  problem  may  be  stated  as  follows: 

These  are  available  N  samples 


vi  •  VJ  vn 

which,  it  is  assumed,  are  known  to  have  come  from  either  the  distribution 


dPl  ■  ve  2  dv 


or  the  distribution 


dPt  •  ve  2  Ig(av)dv 


(210) 


(211) 


the  former  being  the  distribution  of  the  envelope  of  noise  alone,  and  the  latter 
the  distribution  of  signal-plus-noise. 

The  probability  that  all  of  the  variates  ^  •***%  came  from  the  second  dis¬ 
tribution  is  simply 


dPN2  "  dPi{Vl)dP*(Vi)  *“*  dPM 


(212) 
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whereas  the  probability  that  they  all  came  from  the  first  distribution  is 


dPm  m  dPi(\)dPi(vi)  — "  dPi(vs) 


(213) 


The  ratio  of  dPy^  to  dPN[  is  the  best  measure  of  the  likelihood  that  all  the  variants 
came  from  the  signal-plus-noise  distribution.  It  can  be  shown  that  any  monotonic 
function  of  this  ratio  gives  an  equally  good  significance  test.  One  arbitrarily 
picks  a  constant  which  the  ratio  must  exceed  to  say  that  it  shows  that  the  variants 
came  from  the  signal-plus-noise  distribution.  This  constant  determines  the  false 
alarm  time. 


Taking  the  ratio  of  Eq. (213)  to  Eq.(212)  and  substituting  values  from  Eqs.(210) 
and  (211)  gives 


dP, 


TT  "i* 


N1  i«l 


dP, 


,vi 


l.J V 

TTv'T 

i*  I 


or 


(214) 


Jo(av<)  *  (215) 


where  X.  is  the  constant  which  determines  the  false  alarm  time. 


Taking  the  log  of  both  sides  of  Eq.(215)  gives 

j!  log«  Jo(av.)  *  1<J8*  k  *  T  (216> 

i*o 

Note  that  nothing  has  been  said  in  the  foregoing  discussion  about  integration. 
Now,  however,  Eq.(216)  says  that  the  best  thing  to  do  is  take  the  log  of  IQ  of 
each  variate,  add  these  functions  for  each  variate,  and  require  the  sum  to  exceed 
a  certain  value.  Clearly  this  calls  for  a  detector  and  integrator  which  has  the 
combined  law 


y  ■  log  I0(av)  (217) 

The  meaning  of  this  result  is  really  quite  remarkable  (at  least  to  one  who  is  not 
a  statistician).  It  says,  in  effect,  that  by  having  the  sum  only  of  / V  variates 
which  have  been  subjected  to  the  law  y  >  log  I0{av),  one  has  as  much  useful  in¬ 
formation  as  if  the  individual  values  of  each  of  the  variates  were  known  (as  far 
as  determining  to  which  distribution  the  variates  belong)*. 


*  If  the  two  distribution  functions  to  be  distinguished  srs  normal,  then  the  simple  sum  of 
the  N  variates,  or  the  assn,  is  the  best  criterion.  In  other  words,  a  linear  law  would 
be  the  best  if  the  envelopes  of  noise  and  signal-plus-noise  were  norms  1 ly  dis tri buted. 
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Suppose  that  the  signal  strength  is  very  small  (which  would  make  /V  large  for 
any  reasonable  probability  of  detection).  Then  ( av)  2;  1  ♦  aJvJ/4  and 


/  asv3\  a3 v* 
y  ■  log  f0(av)  3s  log  “ — 


(218) 


In  this  case i  the  square  law  is  seen  to  be  the  best  possible  choice.  If,  on  the 
other  hand,  the  signal  strength  is  largo,  I4  ( av)  2;  eav/\/2nav  and 

ettV  1 

y  ■  log  l.(av)  3;  log  ~  .  %  av  -  —  log  27Tav  3s  at>  (219) 

0  v2nav  2 

Thus,  for  large  signals  (usually  small  (V)  the  linear  law  is  best. 

It  should  be  pointed  out  that  the  results  for  the  two  extreme  cases,  square 
and  linear  law,  are  not  very  different  (see  Fig. 42),  and  in  practice  a  linear  de¬ 
tector  would  usually  be  preferred  on  account  of  its  relative  immunity  to  saturation 
by  large  signals. 

In  the  case  of  a  human  operator  it  is  difficult  to  say  what  law  is  used  in 
the  process  of  integration.  Thus  if  a  linear  detector  were  used  in  the  receiver,  it 
is  conceivable  that  the  operator  might  mentally  take  the  sum  of  the  squares  in  his 
integration  process,  with  a  net  over-all  square  law  effect. 


SIGNAL-PLUS-NOISE  MINUS  NOISE  -  LINEAR  LAW 


This  case  is  of  special  interest  because  of  the  method  which  must  be  used  in 
obtaining  the  solution.  Since  the  characteristic  function  for  the  linear  case  cannot 
be  found,  it  is  necessary  to  determine  the  moments  for  a  composite  variate  directly 
from  the  momenta  for  the  signal -plus-noise  distribution  and  those  for  the  noise 
distribution  alone, 


Using  a  double  subscript  notation,  in  which  the  first  index  represents  the 
number  of  the  distribution  function  and  the  second  index  represents  the  order  of 
the  moment,  the  following  formulae  can  be  derived  at  once  by  successive  differen¬ 
tiations  of  the  product  of  the  characteristic  functions  of  the  individual  distribution 
functions: 


h  •  vu  *  vi» 


*2  "  V12  +  U22  +  2VUV22 


^  "  V23  +  V13  *  *{VllV22  *  V!2V2t) 


V4  ■  *  *>4  *  ♦(«'UV33  *  VUV\l)  *  tV22VZ 


(220a) 

(220b) 

(220c) 

(220d) 


V  it  V  ♦  V 


VlSVu)  *  “<*'«*'*  *  V22Vu)  *  i0VliV23 


6(V2SVll 


(220e) 
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The  first  set  of  moments  ere  those  for  one  variate  of  signal  plus  noise  given 
by  Eqs.(109)  and  (llOa-d),  The  second  set  of  moments  are  for  one  negative  variate 
of  noise  alone.  These  are  simply  obtained  from  the  first  set  of  moments  by  putting 
*  ■  0  and  multiplying  the  odd  moments  by  -1. 

The  details  will  not  be  given,  since  the  results  bear  the  same  relation  in 
general  to  the  square  law  case  as  they  do  when  a  noise  variate  is  not  subtracted 
from  each  signal-plus-noise  variate. 


USE  OF  SO-CALLED  DETECTION  CRITERIA 


Lawson  and  Uhlenbeck  have  made  use  of  a  quantity  which  is  the  shift  in  average 
value  of  a  distribution  of  signal-plus-noise  from  that  of  noise  alone  divided  by 
the  standard  deviation  of  noise  alone,  which  they  call  the  detection  criterion.  In 
symbolic  form 


k 


Vl S*H  ~VIN 


(221) 


This  quantity  is  also  called  the  deflection  criterion,  and  it  is  implied  that  it 
must  be  of  the  order  of  unity  or  greater  to  have  a  reasonable  probability  of  detection. 


For  the  square  law  detector,  using  the  results  of  Eqs.(81a-b)  and  (85c),  the 
criterion  becomes 


*  ■  xSFi 


(222) 


and  for  the  linear  detector 


k 


0.957  xVN 


(223) 


assuming  x  to  be  small. 

The  object  of  these  criteria  is  to  show  the  variation  in  necessary  signal* 
to-noise  ratio  as  a  function  of  the  number  of  pulses  integrated.  The  results  for 
k  in  Eqs.(222)  and  (223)  may  be  derived  rigorously  from  the  basic  distribution 
equations  if  the  central  limit  theorem  is  assumed  to  hold  and  for  probability  of 
detection  equal  to  0.50. 

However,  it  is  found  from  the  actual  results  presented  inNo.l,  Figs. 1-50, 
that  the  square  root  of  N  law  given  by  the  detection  criteria  is  not  closely  fol¬ 
lowed,  even  for  N  as  large  as  1,000.  If  a  law  of  the  form 
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is  assumed,  the  exponent  B  may  be  obtained  from  the  data  of  Figs. 1-50,  No.l.  The 
results  are  given  in  Figs. 55  and  56,  No.l.  It  is  seen  that  6  goes  from  1.0  at  N  • 
1  to  around  0.75  at  JV  ■  1,000.  As  pointed  out  earlier  (page  28),  B  »  1  for  co¬ 
herent  integration. 

It  has  been  said  that  the  law  seems  to  fit  observed  data  fairly  well.  It 
is  the  belief  of  the  author  that  this  is  a  coincidence  that  arises  from  the  fact 
that  the  losses  due  to  nonlinear  integration  by  cathode  ray  tubes,  and  human  operator 
losses,  tend  to  just  about  equal  the  difference  between  N*  and  N*  so  that  the  iY^ 
law  actually  seems  to  fit  the  observed  data. 

It  is  rather  interesting  to  note  that  if  the  detector  law  is  assumed  to  be  of 
the  form  y  •  u",  the  detection  criterion  turns  out  to  be 


k  - 


A  graph  of  this  function  shows  a  very  broad  maximum  of  1  at  n  ■  2.  Thus  this  is  a 
special  case,  showing  that  for  large  N  the  square  law  is  the  best  of  the  particular 
class  of  functions  vn.  This  is  not  as  general  as  the  proof  on  page  56  which  shows 
that  the  square  law  is  the  best  of  all  possible  functions  for  small  x. 


COLLAPSING  LOSS  -  INTEGRATION  OF  GREATER  NUMBER 
OF  NOISE  VARIATES  THAN  OF  S IGNAL-PLUS -NOISE  VARIATES 


In  many  radar  applications,  an  additional  number  of  noise  variates  are  integrated 
along  with  a  given  number  of  signal-plus-noise  variates.  Such  is  the  case  when 
three-dimensional  data  are  compressed  onto  a  two-dimensional  presentation,  or  with 
a  C  scope  where  range  is  not  shown.  The  loss  so  occasioned  is  called  a  collapsing 
loss(,#).  An  effect  of  the  same  kind  is  caused  if  the  spot  of  a  cathode  ray  tube 
indicator  moves  less  than  it3  diameter  in  a  pulse  length (11).  Again,  if  the  video 
bandwidth  is  narrow  compared  with  the  IF  bandwidth,  the  same  sort  of  thing  happens. 
All  three  effects  are  handled  by  assuming  a  given  collapsing  ratio,  p,  which  is 
defined  by 
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where 


N  *  number  of  signal-plus-noise  variates  integrated 

It  »  number  of  effective  additional  noise  variates  integrated. 


In  the  case  of  loss  caused  by  low  writing  speed  of  the  cathode  ray  beam,  the  effective 
collapsing  ratio  is  given  approximately  by 


where 


d  m  spot  diameter 
s  ■  writing  speed 
r  ■  pulse  length. 


Where  the  loss  is  caused  by  a  video  amplifier,  the  equivalent  defining  equation  is 


P*i£  “  3 

where  1 


’  (228) 


< 


3^  ■  IF  bandwidth  (or  total  combined  RF  and  IF  bandwidth  where  RF 
amplification  is  used) 

Bv  ■  video  bandwidth. 

Mathematically,  the  treatment  necessary  to  take  account  of  H  extra  noise  variates 
is  rather  simple.  It  is  only  necessary  to  multiply  the  characteristic  function  for 
N  signal -plus -noise  variates  by  the  characteristic  function  for  M  noise-alone  variates. 
In  the  square  law  case,  this  results  in 

^  _  e-Mj) 

(p+ iy’M*  (p+i 

It  is  apparent,  by  comparison  with  Eq.(36),  that  the  results  obtained  for  p  »  1 
can  be  used  directly  to  obtain  results  for  any  p. 

Care  must  be  taken  in  obtaining  the  bias  level,  however.  Without  the  M  extra 
noise  variates,  the  relation  n'  •  n/N  is  used  to  find  the  required  signal-to-noise 
ratio,  the  added  noise  variates,  the  number  of  groups  of  pulses  integrated 

may  oi  r.ny  not  remain  the  same.  In  the  case  of  video  mixing,  where  the  output  of 
two  independent  radars  is  superimposed  on  the  same  indicator,  the  number  of  groups 
of  pulses  integrated  is  constant,  which  means  that  n'  is  constant. 

In  the  other  cases  where  the  loss  is  caused  by  narrow  video  amplifiers,  col¬ 
lapsing  of  coordinates,  or  slow  writing  speed,  the  number  of  independent  groups 
of  pulses  integrated  is  reduced  by  the  factor  peff  so  that  n  remains  constant  as  is 
easily  seen  from  the  equations 

n  ■  n' N  (no  loss)  (230) 

n  -  (pn’)(M+N)  -  n\V  (with  loss).  (231) 


e  P*1 
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(229) 


The  collapsing  loss  is  defined  as 


10  log10 


(232) 


where  is  the  required  signal-to-noise  ratio  with  M  extra  noise  variates,  and 
is  the  signal-to-noise  ratio  required  with  no  extra  noise  variates,  such  that 
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the  probability  of  detection  is  the  same  in  both  cases.  This  fixed  probability 
level  will  usually  be  taken  as  0.90. 

The  procedure,  after  finding  ,  is  to  get  the  required  bias  from  either  Fig. 

8  or  9 ,  depending  on  whether  n  or  n'  is  held  constant,  using  pN  as  the  number  of 
variates.  From  the  cumulative  distribution  functions  graphed  in  Figs.  13  to  32,  the 
value  of  *a  is  found  by  multiplying  the  finding  x  for  pN  variates  to  give  P  *  0.90 
and  multiplying  this  value  of  x  by  p.  The  reason  for  multiplying  by  p  is  apparent 
on  referring  to  £q.(229). 

The  results  of  the  calculation  are  shown  in  Figs.  49  to  52  where  L  is  plotted 
as  a  function  of  N  for  P  «  0.90  and  N  ■  10s.  Also  given  are  curves  of  defined  by 
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It  has  cotnno;.ly  been  said  that  8g  should  be  1/2  <8a>  *  <3#> .  This  statement  is  sometimes 
derived  from  the  detection  criterion  given  on  page  58. 

From  Fig.  57  it  is  seen  that  if  n'  is  constant,  6  does  approach  1/2  as  A- to. 
However,  Q  is  much  smaller  for  reasonably  small  N.  In  the  case  of  n  constant,  the 
square  root  law  is  not  even  approached  as  an  asymptote. 

It  was  found  that  the  values  of  L  and  6  are  only  slightly  dependent  on  the 
original  values  of  n  and  P. 


ANTENNA  BEAM  SHAPE  LOSS 


It  has  so  far  been  assumed  that  the  antenna  pattern  was  flat  over  the  half-power 
beamwidth  and  zero  elsewhere.  In  any  practical  case  the  beam  shape  may  usually  be 
approximated  by  a  Gaussian  curve  which  will  hold  fairly  well  out  to  ±  the  beamwidth 
from  the  point  of  maximum  gain.  In  the  case  of  a  searchlighting  antenna,  the  re¬ 
turned  pulses  will  all  fall  at  the  same  place  in  the  beam,  and  if  this  does  not 
happen  to  fall  at  the  maximum  of  the  beam,  the  loss  may  easily  be  taken  into  account 
by  modifying  the  expression  for  gain  used  in  Eq.(9),  No.l  for  calculating  such  that 


where 


(234) 


9a  »  azimuth  angle  between  target  and  antenna  axis 
9f  ■  elevation  angle  between  target  and  antenna  axis 
Ba  ■  half-power  azimuth  beamwidth 
Bf  «  half-power  elevation  beamwidth 


If  the  antenna  is  scanning,  the  problem  is  entirely  changed  because  the  suc¬ 
cessive  returned  pulses  will  be  of  different  magnitude.  It  is  obvious  that  as  the 
antenna  scans  past  a  target,  pulses  should  be  integrated  out  to  some  point  where 
the  principle  of  diminishing  returns  sets  in.  It  is  not  too  difficult  to  determine 
this  point  and  to  calculate  the  loss  occasioned  due  to  the  beam  shape  as  compared 
with  the  ideal  case(4a).  A  complete  treatment  which  covers  the  general  case  of  delay 
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of  the  received  pulse  relative  to  the  transmitted  pulse,  off  axis  in  elevation  ' 

while  scanning  in  azimuth,  and  random  orientation  of  the  pulse  pattern  relative 
to  the  antenna  pattern  is  quite  involved.  However,  the  solution  of  some  special  ! 

cases  hss  shown  the  general  character  to  be  expected  of  the  results. 

The  integration  of  pulses  should  be  carried  to  about  1. 1  times  the  half-power 
beamwidth.  This  figure  is  practically  independent  of  the  signal  strength  (range) 
and  the  number  of  pulses  per  half-power  beamwidth.  When  the  optimum  number  of  pulses 
are  integrated  there  will  be  an  average  losa  over  the  ideal  case  which  assumes 
constant  gain  between  the  half-power  points.  This  loss  is  in  the  neighborhood  of 
1.5  db  and  does  not  depend  much  on  signal  strength  or  number  of  pulses  per  half¬ 
power  beamwidth.  Since  this  loss  is  so  small  it  was  not  considered  worth  while  to 
reproduce  all  the  detailed  calculations  here. 

It  should  be  mentioned  that  special  care  is  necessary  when  one  considers  rates 
of  antenna  scanning  so  fast  that  about  only  1  hit  per  beamwidth  is  obtained.  In 
this  case  it  may  be  expedient  to  make  the  receiving  antenna  lag  the  transmitting 
antenna  to  compensate  for  the  time  of  travel  of  the  pulse,  or  to  step-scan,  that 
is,  move  the  antenna  in  discrete  steps  rather  than  continuously. 

In  order  to  calculate  the  probability  of  detection  in  any  caae  where  the  sue-  , 

cessive  returned  pulses  have  different  signal  strengths,  it  is  necessary  to  obtain 
the  over-all  characteristic  function  by  multiplying  the  characteristic  functions 
for  each  pulse.  Using  this  method  it  is  not  difficult  to  work  out  the  needed  results 
in  any  particular  case, 

LIMITING  LOSS  ’ 

If  limiting  occurs  anywhere  in  the  receiver,  the  probability  of  detection 
will  be  lowered,  everything  else  being  held  constant.  The  video  amplifier  is  the 
first  place  where  limiting  will  probably  occur.  Let  the  limiting  ratio  be  defined 
as  the  ratio  of  the  limit  level  to  the  R.M.S.  noise  level.  Limiting  can  then  be 
represented  mathematically  by  replacing  the  probability  density  function  at  the 
detector  output  by  an  equivalent  function  below  the  limit  level,  and  a  delta  func¬ 
tion  at  the  limit  level  having  an  area  equal  to  all  of  the  area  of  the  original 
function  to  the  right  of  the  limit  level.  The  moments  can  be  calculated  for  these 
new  functions  (noise  alone  and  signal-plus-noise),  and  the  probability  of  detection  I 

found  by  use  of  the  Gram-Gharlier  series  as  usual.  The  calculations  are  quite  tedious 
and  will  not  be  reproduced  here.  The  main  conclusions  are  that  if  the  number  of 
pulses  integrated  is  large,  the  limiting  loss  is  only  a  fraction  of  a  db  if  the  i 

limiting  ratio  is  as  large  as  2  or  3,  but  if  only  one  or  two  pulses  are  integrated  j 

the  limiting  ratio  must  be  in  the  neighborhood  of  10  to  prevent  a  serious  loss.  ; 

Limiting  in  the  output  of  the  integrator  can  also  cause  a  loss,  but  this  loss  ) 

is  small  compared  to  the  loss  caused  by  limiting  of  the  individual  pulses  in  most 

practical  cases. 


EFFECT  OF  SIGNAL  INJECTION  ON  PROBABILITY  OF  DETECTION 

It  has  been  proposed  that  the  minimum  detectable  signal  can  be  decreased  by 
the  injection  of  an  RF  or  IF  carrier  voltage  that  adds  linearly  to  the  received 
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echo  and  the  receiver  noise <41>.  The  theory  is  that  the  total  signal  will  then  be 
large  compared  with  the  noise,  and  thus  the  so-called  modulation  suppression  that 
occurs  in  the  process  of  detection  with  small  signals  will  be  eliminated. 

In  such  a  process,  the  coherence  of  the  injected  signal  with  the  received 
echo  must  be  taken  into  account.  If  the  target  is  moying,  then  the  successive  re¬ 
ceived  pulses  may  be  considered  to  be  random  in  phase,  so  that  the  injected  signal 
will  necessarily  be  noncoherent  with  the  echo.  Analysis  has  shown  that  in  this  case 
the  probability  of  detection  decreases  continuously  as  the  magnitude  of  the  injected 
signal  increases,  assuming  a  linear  or  square  law  detector.  However,  it  can  be 
shown  that  the  best  possible  detector  law  starts  to  change  radically  as  soon  as 
the  injected  signal  strength  becomes  comparable  to  noise.  The  analysis  of  proba¬ 
bility  of  detection  when  the  detector  function  is  altered  to  take  into  account 
the  injected  signal  has  not  been  completed.  Preliminary  estimates  indicate  that 
there  will  be  only  a  small  decrease  in  sensitivity  in  this  case. 

It  might  be  imagined  that  coherence  could  be  obtained  in  a  system  using  only 
one  hit  per  target  but  having,  say,  20  separate  receiver  channels  with  20  separate 
injection  oscillators  having  phases  spaced  12  degrees  apart..  Thus,  the  return  echo 
would  be  nearly  coherent  with  some  one  of  the  channels.  Theoretically,  the  improvement 
in  this  channel  would  be  about  1  db.  However,  even  this  improvement  would  be  just 
offset  by  the  increased  false  alarm  number  due  to  the  multiple  channels,  so  that 
the  over-all  system  improvement  would  be  nil.  It  seems  that  there  is  no  way  to 
increase  system  sensitivity  to  moving  targets  by  signal  injection. 

Ihere  is  some  possibility  of  increasing  sensitivity  for  stationary  targets 
by  coherent  signal  injection,  but  it  is  difficult  to  imagine  a  practical  situation 
where  such  a  method  would  be  of  any  use. 

PROBABILITY  OF  DETECTION  WITH  MOVING  TARGET  INDICATION  SYSTEMS 

The  analysis  of  the  probability  of  detection  for  MTI  systems  is  quite  compli¬ 
cated.  It  depends  on  the  type  of  receiver  (lin-log  limiting  or  IAGC) ,  the  type 
of  detector,  and  the  characteristics  of  the  storage  device  used.  For  a  nonfluc¬ 
tuating  clutter  and  no  scanning  noise,  the  effect  of  the  clutter  with  or  without 
the  addition  of  a  coherent  oscillator  is  much  the  same  as  that  of  the  injected 
carrier  discussed  in  the  previous  section.  If  a  suitable  detection  system  is  used,  the 
sensitivity  may  be  reduced  by  a  small  amount,  due  to  the  addition  of  the  coho,  perhaps 
by  1  to  3  db. 

The  sensitivity  of  an  MTI  system  for  high  probabilities  of  detection  is  further 
reduced  due  to  the  fact  that  the  target  may  be  moving  at  a  speed  differing  from 
one  of  the  so-called  optimum  speeds.  This  effect  is  quite  complicated  and  is  similar 
to  that  caused  by  a  random  variation  of  the  cross  section  of  a  target  with  aspect. 
A  method  of  quantitatively  treating  these  problems  has  been  developed  and  will  be 
presented  in  detail  in  a  future  report. 

If  there  is  a  rluctuation  component  in  the  clutter,  due  either  to  the  movement 
of  the  clutter  itself  or  to  the  scanning  of  the  antenna,  the  effect  will  be  to 
increase  the  amount  of  noise  at  the  receiver  input.  This  can  be  taken  into  account 
by  an  appropriate  adjustment  in  the  value  of  the  noise  figure  of  the  receiver  that 
will  change  by  the  correct  amount. 
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TABLES  OF  THE  DERIVATIVES  OF  THE  ERROR  FUNCTION 


In  order  to  make  efficient  use  of  Gram-Charlier  series,  it  is  necessary  to 
have  a  good  table  of  the  derivatives  of  the  error  integral  (the  <t>  functions  of  Eq. 
62).  No  satisfactory  table  was  in  existence  at  the  time  this  report  was  written. 
Typical  of  the  available  tables <a>  were 

Fry<*>  n  ■  1(1)6,  *  ■  0(.1)4  5  decimals 

Jorgensen  n  ■  1(1)6,  *  ■  0(.01)4  7  decimals 

and  an  unpublished  table  of  the  W.P.A. ,  giving 

n  “  1(1)14,  *  ■  0(. 1)8.4  20  decimals 

RAND  therefore  decided  to  calculate  a  suitable  table  with  the  aid  of  its  IBM 

equipment.  This  has  resulted  in  a  table  of  Hermite  polynomials,  as  well  as  in  the 

derivatives  of  the  error  integral,  giving 

n  «  1(1)10,  *  ■  0(. 01)12.0  6  significant  figures 

A  limited  number  of  these  tables  are  available  at  the  present  time.  ( RAND  Document 
D-350,  A  Table  of  Hermite  Polynomials  and  the  Der  ivatives  of  the  Error  Function. ) 
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